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ABSTRACT 


A  critical  reexamination  is  made  of  existing  experimental  data  on 
the  velocity,  temperature,  and  concentration  profiles  in  a  low-speed, 
constant  property  turbulent  boundary  layer  with  surface  mass  injection. 

It  is  found  that  a  model  of  the  boundary  layer  can  be  constructed  that  fits 
the  local  profile  data  as  well  as  the  local  drag,  mass,  and  heat-transfer 
coefficients  at  the  surface.  The  essential  features  of  the  model  include 
a  "law  of  ..*ie  wall"  that  provides  a  continuous  and  smooth  profile  from 
the  surface  to  the  fully  turbulent  region  and  a  "wake  region"  that  joins 
the  inner  region  smoothly  and  is  characterized  by  an  eddy  diffusivity  that 
is  constant  with  respect  to  distance  from  the  wall.  Four  empirical  param¬ 
eters  appear  in  the  analysis:  the  Prandtl  mixing  length  constant,  K;  the 
Clauser  eddy  diffusivity  Reynolds  number;  an  inner  region  mixing  length 
constant  analogous  to  Rotta's  intercept  constant;  and  a  value  of  turbulent 
Schmidt  or  Prandtl  number.  Constant  values  of  these  parameters  are  found 
to  be  suitable  in  predicting  data  over  a  wide  range  of  surface  blowing  or 
suction  and  Reynolds  number. 
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MASS,  MOMENTUM,  AND  HEAT  TRANSFER  WITHIN  A  TURBULENT 
BOUNDARY  LAYER  WITH  FOREIGN  GAS  MASS  TRANS FFR 
AT  THE  SURFACE 


1 .  INTRODUCTION 

The  most  effective  current  techniques  for  thermally  protectinq  missile 
nose  cones  or  manned  reentry  vehicles  from  aerodynamic  heating,  namely,  tech¬ 
niques  involving  ablation  or  transpiration  cooling,  are  characterized  by  a 
mass  transfer  of  gas  through  the  surface  of  the  body.  Because  of  this  mass 
transfer  and  the  attendant  surface  roughness  of  ablating  surfaces  or  porous 
skins,  turbulent  boundary  layers  can  be  expected  to  exist  on  the  surfaces  of 
reentry  vehicles  even  at  the  Reynolds  numbers  corresponding  to  rather  shallow 
reentry  angles.  Thus,  an  understanding  of  the  behavior  of  turbulent  boundary 
layers  under  the  conditions  of  reentry  and  with  simultaneous  surface  mass 
transfer  is  of  considerable  importance  in  the  design  of  reentry  body  thermal 
protection  systems.  In  addition,  this  understanding  is  essential  in  the  pre¬ 
diction  of  electromagnetic  radiation  emission  or  of  radar  reflectivity  from 
the  boundary  layer  or  the  initial  stages  of  the  wake. 

When  an  examination  is  made  of  the  published  literature  on  this  subject, 
it  becomes  apparent  that  rather  tenuous  extrapolations  of  tests  of  limited 
scope  and  of  early  heuristic  theories  have  been  made  to  predict  the  behavior 
of  the  dynamically  complex,  chemically  reacting  flows  about  the  reentry  bodies. 
The  work  contained  in  this  report  constitutes  a  reexamination  of  the  problem, 
with  the  intent  of  arriving  at  a  model  of  the  t  rbulent  boundary  layer  that 
is  entirely  consistent  with  data  obtained  under  successively  more  complex 
situations . 

l.i  Historical  Review 

Before  a  review  is  made  of  the  historical  background  of  the  current 
research,  it  will  be  worthwhile  to  present  a  description  of  the  general  struc¬ 
ture  of  the  turbulent  boundary  layer  as  it  is  presently  understood.  This 
will  provide  the  reader  with  the  hindsight  the  authors  enjoyed  while  perform¬ 
ing  this  c**  tical  reexamination,  and  will  clarify  some  of  the  approaches 
taken  in  the  analysis.  The  structure  of  the  turbulent  boundary  layer  is  de¬ 
scribed  in  detail  in  References  1  and  2  and  is  summarized  in  the  following 
sections . 


1-2 


1.1.1  General  description  of  the  turbulent  boundary-layer  structure 

Generally,  the  boundary  layer  can  be  considered  to  be  diviied  into  four 
successive  layers  from  the  surface  to  the  outer  edge. 

(a)  The  layer  closest  to  the  surface  is  the  "laminar  sublayer" 

or  the  "viscous  layer."  This  layer  is  characterized  by  shear  stresses  that 
arc  predominately  caused  by  molecular  viscosity.  The  flow  here  need  not  be 
steady,  but  can  fluctuate  as  a  result  of  turbulence  in  the  layers  above. 

Where  no  mass  transfer  occurs  normal  to  the  surface,  the  time-averaged  shear 
in  this  layer  is  essentially  constant  with  distance  from  the  surface. 

(b)  The  next  layer  away  from  the  surface  is  the  "buffer  layer." 

In  this  layer,  the  shear  stress  is  dependent  on  both  the  molecular  viscosity 
and  the  turbulent  motion  (Reynolds  stresses).  Again,  with  no  surface  mass 
transfer,  the  time-averaged  shear  stress  in  this  layer  is  essentially  con¬ 
stant  with  distance  from  the  surface.  The  term  "buffer  layer"  was  employed 
originally  by  von  Karman  to  indicate  the  presence  of  this  layer  between  a 
completely  laminar  and  a  completely  turbulent  layer. 

(c)  The  third  layer  from  the  surface  is  called  "fully  turbulent  region. 
In  this  layer,  the  shear  stress  is  dependent  primarily  on  the  turbulent 
motions.  The  velocity  distributions  in  this  layer  have  been  shown  experimen¬ 
tally  to  be  dependent  on  che  local  wall  conditions,  and  rather  insensitive 

to  the  local  boundary-layer  thickness  or  pressure  gradients.  This  has  led 
Coles  to  apply  the  term  "law  of  the  wall"  to  describe  this  layer,  together 
with  layers  (a)  and  (b) ,  above,  References  3  and  4. 

(d)  The  fourth  layer  is  also  fully  turbulent;  however,  this  layer's 
characteristics  are  controlled  more  by  the  boundary-layer  thickness  than  the 
local  shear  at  the  wall.  The  flow  in  this  region  is  generally  characterized 
by  eddies  that  are  similar  to  those  that  exist  in  wakes  and,  thus,  Coles 
could  represent  the  velocity  in  this  region  for  nonporous  walls  by  a  "law  of 
the  wake,"  Reference  4.  Others,  References  5  and  6,  have  applied  the  term 
"velocity  decrement  law"  to  velocity  correlations  for  this  same  region. 

The  relative  thickness  of  the  four  layers  described  above  are  indicated 
in  Figure  1,  reproduced  from  Reference  2.  The  values  shown  in  this  figure 
correspond  to  a  turbulent  boundary  layer  on  a  nonporous  flat  plate  in  low- 
speed  flow,  where  the  properties  in  the  boundary  layer  are  essentially  con¬ 
stant,  and  where  the  Blasius  skin-friction  law  is  applicable.  It  reveals 
that  the  relativt  thicknesses  of  the  inner  layers  diminish  with  increases  in 
the  local  length  Reynolds  number.  At  Re^  •*  108 ,  the  boundary-layer  law  of 
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the  wall  applies  over  only  5  percent  of  the  thickness  of  the  boundary  layer, 
while  the  wake  region  occupies  95  percent  of  the  boundary  layer.  When  con¬ 
sideration  is  made  that  heat  and  mass  transfer  within  the  boundary  layer 
with  Prandtl  numbers  and  Schmidt  numbers  less  than  unity  result  in  larger 
thicknesses  of  the  diffusional  boundary  layers  relative  to  the  flow  or 
momentum  exchange  boundary  layers  ,  the  presence  of  the  wake -type  flow 
becomes  even  nore  important  and  requires  consideration  at  lower  Reynolds 
numbers.  Also,  surface  injection  tends  to  increase  the  boundary-layer  thick¬ 
ness  which  may  force  a  consideration  of  the  wake  layers  at  lower  length 
Reynolds  numbers  than  indicated  in  Figure  1. 

1.1.2  Early  contributions 

In  the  late  1940's  and  early  1950's,  a  need  arose  for  more  effective 
cooling  systems  for  missile  surfaces  experiencing  aerodynamic  heating  through 
turbulent  boundary  layers.  It  was  known  then,  from  research  with  laminar 
boundary  layers,  that  "transpiration"  of  gas  at  a  surface  was  a  very  effec¬ 
tive  cooling  mechanism,  Reference  7.  Because  a  qualitatively  similar  effec¬ 
tiveness  was  expected  to  occur  with  the  turbulent  boundary  layer,  research 
began  on  this  subject  at  MIT,  and  then  independently,  a  few  years  later, 
at  Convair  and  the  Ames  Lat oratory  of  the  NACA  (now  NASA). 

The  MIT  work,  under  the  direction  of  H.  S.  Mick ley,  was  largely  experi¬ 
mental  during  its  early  periods,  Reference  8.  This  excellent  work  has  con¬ 
tinued  to  this  day,  References  9  through  15  and  17  through  20,  and  provides 
the  bases  for  many  of  the  deductions  of  this  report  that  are  discussed  in 
detail  later.  Reference  9  served  as  the  basis  for  the  empirical  parameters 
of  the  early  heuristic  theories  described  in  the  following  paragraphs. 

To  satisfy  the  needs  for  effective  missile  cooling  s_  -jms ,  the  specula¬ 
tive  analyses  of  References  21  and  22  were  performed  to  yield  indications  of 
the  effect  of  surface  mass  transfer  of  air  under  supersonic  conditions  where 
the  tnen  existent  MIT  data  of  Reference  9  did  not  apply  directly.  Reference 
21  employed  the  classical  Prandtl  mixing-length  theory,  with  variable  density, 
across  the  entire  boundary  layer  on  a  flat  plate.  Thus,  the  viscous  sub¬ 
layer,  the  buffer  layer,  and  the  wake  regions  were  not  considered-  In  addi¬ 
tion,  both  the  laminar  and  tu’  °nt  Prandtl  numbers  were  set  equal  to  unity 
in  the  energy  equation.  The  local  shear  stress  distribution  in  the  boundary 
layer  was  modified  to  account  for  the  surface  mass  transfer;  that  is,  the 
local  shear  stress  was  expressed  as 
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The  omission  of  the  wake  region  resulted  in  the  mathenat ical  requirement  of 
a  continual  increase  of  shear  stress  with  distance  away  from  the  wall1  even 
with  the  knovn  boundary  condition  that  the  average  shear  stress  must  approach 
zero  at  the  cuter  edge  of  the  boundary  layer.  Implicit  here  was  the  hope 
that  the  equally  erroneous  ever  increasing  mixing  length,  1  =  ky,  somehow 
compensated  to  yield  adequate  velocity  distribution  neai  the  outer  edge  cf 
the  boundary  layer.  The  closed-form  solution  for  the  velocity  distribution 
in  the  turbulent  region  obtained  in  this  analysis  could,  by  an  adjustment  of 
the  arbitrary  constants  of  integration,  be  made  to  fit  the  compressible  analog 
of  the  classic  logarithmic  velocity  distribution  for  the  case  of  zero  blowing, 
that  is,  Reference  23.  Thus,  an  additional  implied  assumption  was  made  that 
the  "constants "  of  integration  were  not  dependent,  themselves,  on  the  blowing 
rate . 

An  analysis  based  on  the  identical  physical  assumptions  as  in  Reference 
21  was  perforned  independently  in  Reference  24,  except  that  the  latter  inves¬ 
tigation  was  restricted  to  constant  fluid  properties.  Again,  Reference  9 
was  used  as  the  experimental  basis  for  comparison.  It  was  found  that  the 
theory  generally  predicted  the  trends  in  the  data  on  boundary-layer  growth 
and  skin  friction  (to  ~  25  percent) ;  however,  the  predicted  and  measured 
velocity  profiles  disagreed  considerably  in  shape. 

To  account  for  the  Prandtl  Timbers  other  than  unity  in  the  heat-transfer 
rates,  Reference  22  added  a  laminar  sublayer  to  the  fully  turbulent  outer 
region,  again  characterized  by  the  . randti  mixing-length  theory  and  the  shear 
distribution  represented  in  Equation  (1-1).  The  omission  of  a  buffer  layer 
caused  a  discontinuity  in  the  slope  of  the  velocity  profile  to  occur  at  the 
interface  of  tie  sublayer  and  fully  turbulent  regions.  It  was  found  that  the 
thickness  of  the  sublayer  was  an  important  parameter,  and  recourse  to  the 
l}w-speed  data  of  Reference  9  was  made  to  define  a  realistic  hypothesis  for 
the  sublayer  thickness.  The  results  of  Reference  22,  while  yielding  quanti¬ 
tative  results  on  the  expected  reduction  in  surface  heat-transfer  rates  that 
agreed  fairly  well  with  subsequent  experimental  data,  Reference  25,  never 
yielded  satisfying  results  concerning  the  velocity  and  temperature  profiles 
over  the  entire  boundary  layer  or  of  the  tempeiature  or  enthalpy  recovery 
factor.  The  theory  showed  an  ever  increasing  value  of  recovery  factor  with 

The  term  u  in  Equation  (1-1)  increases  monotonically  with  distance  from 
the  surface. 
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increasing  length  Reynolds  number  that  was  not  evident  in  the  zero  and  finite 
lowing-rate  data.  Significantly,  no  consistent  Mach  number  effect  was  evi¬ 
dent  in  the  theoretical  results,  and  the  tests  of  Reference  25  did  not  con¬ 
tribute  information  regarding  these  Mach  number  effects,  as  they  were  con¬ 
ducted  at  a  single  Mach  number,  M  -  2.7.  Later  it  was  indicated,  Reference  26, 
that  the  theory  of  Reference  22  agreed  to  about  +25  percent  with  the  skin- 
friction  d. ta  obtained  at  M  =  5.1  from  boundary-layer  velocity  surveys. 

This  agreement  appears,  i.i  retrospect,  to  have  been  fortuitous  because  the 
skin-friction  data  were  obtained  from  an  extrapolation  of  velocity  distribu¬ 
tion  measurements  toward  the  surface,  usually  a  rather  inaccurate  technique 
with  probes  as  large  as  the  conventional  impact  pressure  tubes  (see  Analysis). 

Theories  similar  to  Reference  22  were  developed  for  the  case  of  injec¬ 
tion  of  a  single  foreign  gas  into  undissociated  air,  References  27  and  28. 

These  theories  were  further  restricted,  however,  to  changes  in  the  physical 
properties  of  the  binary  mixture  caused  only  by  concentration  variations  of 
the  individual  gases  and  did  not  include  chemical  reactions,  temperature,  or 
compressibility  effects.  In  general,  these  theories  indicated  the  marked 
effectiveness  of  light  gas  to  reduce  heat  transfer  end  skin  friction  at  tne 
surface  but,  in  general,  possessed  the  shortcomings  of  Reference  22.  Further, 
Reference  28  only  yielded  the  relationship  between  heat  transfer  and  skin 
friction,  but  not  an  explicit  value  of  either  term. 

When  the  skin-friction  data  of  Reference  29  were  obtained,  it  was  found 
that  the  results  agreed  qualitatively  with  the  theory  of  Reference  27,  in 
that  the  transpiration  of  lighter  gases  was  more  effective  in  reducing  skin 
friction.  The  data  also  showed  that  the  theories  of  References  21  and  22 
failed  to  predict  quantitatively  the  Mach  number  effect  indicated  in  the  data, 
and  that  the  agreement  of  the  predictions  of  Reference  22  with  earlier  data 
was  fortuitous,  as  the  oata  were  obtained  at  a  Mach  number  of  2.7  where 
agreement  happened  to  occur,  see  Figure  2  reproduced  from  Reference  29.  This 
latter  effect  occurred  because  the  data  of  Reference  9,  the  basis  of  arbitrary 
constants  in  Reference  22,  showed  a  lesser  reduction  in  skin  friction  with 
gas  injection  than  the  later  data  of  Reference  15.  The  data  of  Reference  26, 
therefore,  are  inconsistent  with  these  general  results. 

An  analysis  by  Ness,  Reference  31,  extended  the  work  of  Reference  27  to 
include  variations  of  the  properties  in  the  boundary  layer  caused  by  tempera¬ 
ture  variations  as  well  as  by  concentration  distributions.  The  gases  con¬ 
sidered  were  undissociated  air  and  a  light,  chemically  inert  injected  foreign 
gas.  Ness  used  the  two-layer  formulation  of  References  22  and  27  with 
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essentially  the  same  sublayer  thickness  criterion.  Numerical  values  of  sur¬ 
face  skin  friction  and  heat-transfer  rates  were  evaluated  for  the  particular 
case  of  helium  injection  into  air  up  to  Mach  numbers  of  4. 

For  light  gas  injection  where  Schmidt  numbers  can  vary  between  about 
0.2  to  1.6  in  the  sublayer,  Ness  argues  that  the  turbulent  Schmidt  number 
must  be  equal  to  the  Schmidt  number  based  on  molecular  transport  properties 
at  the  interface  of  the  two  layers,  Reference  32.  His  argument  is  based  on 
the  mathematical  result  that  if  the  turbulent  Schmidt  number  is  a  fixed  value 
lying  between  the  extreme  values  of  the  Schmidt  number  in  the  sublayer,  a 
maximum  in  the  concentration  of  the  light  gas  at  the  wall  occurs  when  the 
velocity  ratio  at  the  edge  of  the  sublayer  is  varied  at  constant,  but  large, 
values  of  the  injection  parameter  C,  .  r~  The  concentration  at  the  wr-ll  varies 
monotonically  if  Ness'  criterion  cited  above  is  applied.  Why  Ness  finds  a 
maximum  in  concentration  of  the  light  gas  at  the  wall  with  varying  sublayer 
velocity  ratio  an  unacceptable  anomaly  is  not  -rlear.  As  the  velocity  ratio 
is  increased  with  increased  skin  friction,  the  resistance  of  the  boundary 
layer  to  mass  transfer  tends  to  increase,  causing  a  rise  in  surface  concen¬ 
tration  for  the  increased  flow  rates  corresponding  to  constant  f, .  On  the 
other  hand,  the  increased  flow  rate  also  results  in  an  increased  light  gas 
concentration  at  the  edge  of  the  sublayer  which  increases  the  molecular 
Schmidt  number  in  the  sublayer  that,  in  turn,  tends  to  drop  the  boundary- 
layer  resistance  to  mass  transfer.  Thus,  two  simultaneous,  but  opposite, 
effects  are  occurring  so  that  it  is  not  inconceivable  chat  a  maximum  in  con¬ 
centration  actually  occurs.  Thus,  Ness'  argument  cannot  be  accepted  a  priori 
but  must  be  tested  against  the  concentration  profile  data  that  exist  It 
will  be  shown  later  that  the  assumption  of  a  laminar  Schmidt  number  equal 
to  0.2  and  a  turbulent  Schmidt  number  equal  to  0.75  yields  excellent  agree¬ 
ment  with  concentration  profile  data;  and  since  this  is  in  contradiction  to 
Ness'  criterion  for  the  value  of  the  turbulent  Schmidt  number,  it  can  be 
concluded  that  his  criterion  is  not  a  necessary  condition. 

In  view  of  the  above  discussion,  some  caution  must  be  exercised  in 
accepting  Ness'  numerical  results  which  indicate  a  correspondence  to  the  Mach 
number  effect  found  in  Reference  29.  In  fact,  Ness'  analysis  reduces  tr  that 
of  Reference  22  for  air  injection,  where  little,  if  any,  Mach  number  effect 
is  noted,  see  Figure  2.  It  is  believed  that  the  results  of  Reference  31  may 
be  fortuitous  and  introduced  through  the  arbitrary  choice  of  the  variation  of 

"The  injection  parameter 

J  r  u  u 

C  "  ^  uJc°/2) 


v 


1-7 


i > 


r 

i. 


I 

i 

I  i 

i 

j 


the  turbulent  Schmidt  number.  Again,  Reference  31  possesses  all  the  defects 
cited  earlier  for  two-layer  models  of  the  turbulent  boundary  layer;  for 
example,  ever  increasing  temperature  recovery  factor  with  increasing  Reynolds 
number. 

Regardless  of  the  limitations  of  the  data  cited  and  of  the  two-layer 
theories,  the  practical  r.^eds  in  the  design  of  ablation  bodies  have  led  to 
the  use  of  results  of  these  investigations  to  define  the  heat  flow  blocking 
effect  by  the  injected  gases,  for  example,  Reference  33. 

1.1.3  Recent  contributions 

The  recent  publication  by  Pappas  and  Okuno,  Reference  34,  extends  the 
work  cf  Reference  29  to  local  heat-transfer  rate  measurements  while  helium, 
air,  and  Freon-12  are  injected  from  the  surface  of  a  15°  included  angle  cone. 

The  results  of  this  test  indicate  further  anomalies  not  evident  in  the  pre¬ 
viously  discussed  theories.  These  can  be  summarized  as  follows: 

(a)  No  particular  Mach  number  effect  is  shown  in  plots  of  St/StQ 
versus  F0v0/p  ux  St  for  air  injection.  This  agrees  with  the  theory  of 
Reference  27,  yet  the  skin-friction  data  of  Reference  29  do  not  agree  with 
the  same  theory.  The  implication  here  if  that  the  theoretical  prediction  of 
the  modified  Reynolds  analogy  is  not  correct. 

(b)  The  data  ;%how  little  variation  of  dropoff  in  Stanton  number  v  ith 
Reynolds  number  which  is  inconsistent  with  the  theory  of  Reference  27  for 
helium  injection. 

(c)  At  subsonic  speeds,  the  recovery  factor  for  helium  injection  increases 
from  the  zero  injec  ion  value  of  0.84  to  a  maximum  value  near  1.35  ano  then 
decreases  with  increasing  injection.  At  supersonic  speeds,  the  recovery  fac¬ 
tors  for  all  the  injection  gases  initially  decrease  from  the  zero  injection 
value  and  then  increase  at  higher  injection  rates  to  va'ues  greater  than  unity. 
None  of  these  effects  are  evident  in  the  existent  theories.  It  i s  believed 
(Ref.  35)  that  thermal  diffusion  in  the  sublayer  may  be  causing  these  recovery- 
factor  effects. 

The  general  result  of  these  data  is  that:  the  light  gases  are  not  relatively 
as  effective  in  reducing  heat-transfer  rates  as  has  hitherto  been  expected 
from  the  theories  cited  above.  These  data  are  not  entirely  consistent  with 
the  heat-transfer  measurements  of  Bartle  and  Leadon,  Refeience  36,  where  a 
Mach  number  effect  on  heat-transfer  data  is  still  noted.  For  injected  gases 
other  than  nitrogen,  there  is  insufficient  detail  in  Reference  36  to  make 
direct  comparison  with  the  data  of  Pappas  and  Okuno. 
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In  general,  it  can  be  deduced  from  the  survey  cited  above  that  several 
major  discrepancies  exist  between  the  data  obtained  to  date  and  existing 
theory  dealing  with  air  or  foreign  gas  injection  into  a  turbulent  air  bound¬ 
ary  layer.  In  addition,  certain  unresolved  differences  exist  between  diifer- 
ent  sets  of  data.  The  use  of  these  theories  to  predict  the  behavior  of  tur¬ 
bulent  boundary  layers  on  ablating  reentry  bodies  requires  extrapolation  that 
borders  on  mere  guessing. 

It  is  the  intent  of  the  current  research,  therefore,  to  create  a  newer 
mde’  of  turbulent  boundary  layer  that  is  consistent  with  successively  more 
complex  experimental  results.  In  the  following  paragraphs  the  broader  aspects 
of  this  model  will  be  described  in  order  that  the  details  of  the  analysis 
presented  in  Section  2  can  be  viewed  in  their  proper  perspective. 

As  has  been  previously  noted  (Section  1.1.1),  the  boundary  layer  can  be 
divided  into  essentially  four  parts,  the  first  three  of  which  are  character¬ 
ized  by  the  "law  of  the  wall."  Because  the  characteristics  of  this  portion 
of  the  boundary  layer  have  been  found  to  be  almost  exclusively  dependent  on 
wall  conditions,  it  is  possible  to  formulate  *-.s  solution  in  terms  of  a  "one¬ 
dimensional"  initial  value  problem.  The  initial  values  are  those  that  char¬ 
acterize  the  conditions  at  the  wall  which,  in  general,  vary  from  one  point 
on  the  wall  to  another.  The  term  "one-dimensional"  is  a  ropriate  in  the 
mathematical  sense  that  the  solution  proceeds  out  (in  _>ne  dimension,  y)  from 
a  specified  wall  point. 

The  fourth  region  of  the  boundary  layer  (as  enumerated  in  Section  1.1.1) 
is  strongly  dependent  on  the  history  of  the  boundary  layer  and  on  local  free- 
stream  conditions  and  gradients.  Thus,  the  solution  of  this  portion  of  the 
boundary-layer  problem  is  characterized  by  a  "two-dimensional"  boundary  value 
problem.  The  boundary  values  are  conditions  at  some  initial  plane  normal  to 
the  flow  direction,  conditions  on  a  plane  parallel  to  the  wall  in  the  free 
stream  beyond  the  extent  of  the  boundary  layer,  and  conditions  at  the  inter¬ 
face  between  this  region  and  the  inner  "law  of  the  wall"  region. 

In  the  analysis  presented  in  the  next  section,  these  two  regimes  will 
be  studied,  first  individually  and  then  in  concert,  to  establish  a  "solution" 
for  the  complete  boundary  layer.  The  establishment  of  this  "solution"  is 
dependent,  in  large  measure,  on  the  physical  assumotions  which  will  be  made 
to  obtain  ~  tractable  set  of  mathematical  relations  for  each  of  these  two 
regimes.  These  aia  briefly  summarized  in  the  following  table: 
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Assumptions 

Mixing  length  representative  of 
the  variation  of  turbulent  trans¬ 
port  properties. 

Rate  of  increase  of  mixing  length 
based  on  local  properties  in  the 
boundary  layer 

Molecular  transport  effects  neg- 
1 igible . 

Turbulent  transport  properties 
are  functions  of  the  streamwise 
coordinate,  x,  only. 

Clauser  parameter,  e^'u^*, 
invariant . 

Invariant  turbulent  Schmidt  and 
Prandtl  numbers. 

The  advantages  realized  by  dividing  the  solution  into  two  distinct  regimes 
are  manifold,  if  consideration  is  given  to  the  complexities  of  a  single  unified 
solution.  It  is  indeed  fortunate  that  such  a  division  not  only  reduces  an 
otherwise  intractable  mathematical  problem  into  two  reasonable  straightforward 
problems,  but  also  is  consistent  with  tie  physical  character  of  the  boundary 
layer.  The  assumptions  used  in  the  wake  regime,  although  reasonably  valid 
therein,  are  totally  invalid  in  the  law  of  the  wall  regime.  Likewise  the 
"one-dimensional"  assumption  of  the  law  of  tie  wall  regime,  although  good 
therein,  cannot  be  extended  into  the  wake  regime.  Thus,  the  complexities  of 
one  regime  are  avoided  in  the  solution  of  the  other,  resulting  in  a  profound 
reduction  in  the  overall  complexity  of  the  total  problem. 

In  previous  work,  boundary-layer  studies  have  emphasised  the  "law  of  the 
wall"  regime,  usually  ignoring  the  contributions  of  the  wake  regime.  The 
principal  difference  between  these  analyses  has  involved  the  chosen  represen¬ 
tation  of  the  variation  of  the  turbulent  transport  effects.  In  the  present 
work  these  various  techniques  have  been  evaluated  in  light  of  the  existing 
data  with  a  particular  emphasis  on  their  application  under  conditions  of  sur¬ 
face  mass  addition.  The  choice  of  a  "mixing  length"  model  was  based  primarily 
or.  its  initial  successful  application  under  these  conditions.  The  selected 
analytical  representation  of  the  mixing  length,  although  original,  is  a  con¬ 
sequence  of  a  generalization  and  rationalization  of  existing  representations. 
3ased  on  tne  data  evaluated,  the  two  constants  of  this  representation  do  not 
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vary  with  surface  mass  addition  or  Schmidt  (or  Prandtl)  number  nor  witn  appli¬ 
cation  to  velocity,  temperature,  or  concentration  gradients.  Because  the 
consequent  mathematical  formulation  is  an  initial- value  problem,  its  solution 
is  conceptually  simple. 

Solutions  for  the  outer  porcion  of  the  boundary  layer  (the  wake  regime), 
with  one  notable  exception,  do  not  exist  per  so.  There  are  various  velocity 
profile  generalizations  which  have  been  evolved  and  successfully  applied  (e.g., 
Coles,  Ref.  3;  Black  and  Sarnecki,  Ref.  37;  and  Kendall,  Ref.  17;  the  latter 
two  being  inspired  by  the  former)  but  these  techniques  are  not  supported  by 
any  particular  physical  rationale  relating  to  the  turbulent  transport  proper¬ 
ties  which  dictate  the  profiles.  The  exception  is  the  work  of  Clauser  (Ref.  6). 
Based  on  observations  with  regard  to  the  variation  of  the  turbulent  viscosity 
within  the  boundary  layer,  he  proposed  a  simplification  and  a  generalization, 
namely,  that  the  turbulent  viscosity  coefficient  was  e  function  of  the  stream- 
wise  coordinate  only,  and  that  this  functional  dependence  could  be  specified 
as 


constant 


that  is,  the  turbulent  viscosity  is  proportional  to  the  product  of  the  free- 
stream  velocity  and  the  boundary-layer  displacement  thickness.  In  his  original 
work,  Clauser  restricted  his  solution  to  velocity  profiles  on  impervious  walls. 
In  the  analysis  to  follow,  the  technique  has  been  extended  to  concentration 
and  temperature  profiles  as  well  as  velocity  profiles,  all  under  condition  of 
surface  mass  addition.  The  procedure  involves  the  formulation  of  a  similarity 
solution,  much  akin  to  the  classical  solutions  of  the  laminar  boundary  layer. 
Some  of  the  requisites  for  a  truly  similar  (turbulent)  boundary-layer  solution 
are  found  to  be  more  palatable  than  their  conterparts  in  the  laminar  boundary 
layer.  In  particular,  a  constant  value  f  the  transpiration  ratio  P0V0/P1V1 
is  required  rather  than  its  inverse  decay  with  boundary-layer  thickness  as 
prescribed  for  laminar  similarity  solutions. 

The  remaining  requirements  of  the  solution  of  the  entire  boundary  layer 
involve  the  coupling  of  the  inner  and  outer  solutions  and  generalizations  of 
the  results  obtained  for  specified  boundary  conditions.  This  is  again  an 
almost  purely  mathematical  effort  and  requires  no  further  physical  assumptions, 
beyond  those  already  discussed.  A  total  of  four  invariant  constants  are  util¬ 
ized:  the  Clauser  parameter,  the  turbulent  Prandtl  (or  Schmidt)  nuror'r,  and 

the  two  constants  of  the  mixing  length  expression. 


''■nmifai i 


1-11 


In  checking  for  consistency  between  theory  and  experiment,  emphasis  has 
been  placed  on  velocity,  temperature,  and  gas  specie  concentration  profiles 
in  addition  to  overall  surface  transfer  coefficients,  such  as  the  skin  fric¬ 
tion  coefficient  or  Stanton  number.  The  work  reported  in  this  report  is 
limited  to  the  case  where  the  fluid  properties  '.n  the  boundary  layer  are 
essentially  constant.  The  excellent  agreement  of  the  proposed  boundary-layer 
model  with  much  of  the  data  obtained  at  low  speeds  and  with  small  quantities 
of  foreign  gas  injection  is  quite  encouraging.  Even  the  initial  objectives 
of  this  program  cannot  be  considered  to  be  complete,  however,  until  the 
boundary-layer  model  is  adapted  to  be  consistent  with  existing  data  on  the 
following: 

(a)  Boundary-layer  velocity  profiles  obtained  at  supersonic  speeds  with 
air  injection  at  the  surface. 

(b)  Surface  transfer  coefficients  in  supersonic  flow  with  foreign  gas 
injection  in  moderate  and  large  quantities. 

(c)  Surface  transfer  coefficients  deduced  in  sonic  or  supersonic  flows 
where  intense  chemical  activity  and  high  surface  mass  transfer  takes  place. 

Ultimately,  the  theoretical  model  of  the  boundary  layer  will  have  to  be  tested 
against  future  experiments  in  high-speed  chemically  reacting  boundary  layers 
where  local  profile  data  are  obtained. 
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2.  ANALYSIS 

2.1  Basic  Boundary-Layer  Equations 

Because  of  the  complexity  of  the  ultimate  goals  of  this  study,  it  is 
imperative  that  the  considerably  more  restrictive  objectives  of  the  initial 
work  reported  herein  be  brought  into  proper  perspect ive .  To  aid  in  the 
establishment  of  this  perspective,  this  section  w  11  present  the  basic 
boundary-layer  equations  which  are  appropriate  to  the  general  solutions 
ultimately  sought.  A  sequence  of  restrictions  will  subsequently  be  imposed 
upon  these  general  relations  in  order  to  arrive  at  the  basic  set  of  relations 
which  will  be  considered  in  detail  in  this  report.  In  the  various  develop¬ 
ments  which  will  be  presented  in  the  subsequent  sections  of  this  report, 
means  of  relaxing  many  of  these  restrictions  will  be  discussed.  Initially, 
therefore,  the  boundary- layer  relations  presented  will  include  effects  of 
compressibility,  chemical  reactions  between  multiple  components,  unequal 
binary  diffusion  coefficients,  thermal  diffusion  and  its  reciprocal  partner, 
the  DuFour  effect,  and  temperature  and  concentration-dependent  physical 
properties . 

2.1.1  General  two-dimensional  mass,  momentum,  and  energy  equations 

In  the  derivation  of  the  equations  presented  here,  the  standard  defini¬ 
tions  of  time-averaged  turbulent  quantities  and  relative  order  of  magnitude 
are  employed,  for  example,  References  23  and  38.  Because  the  turbulent 
transport  terms  are  expressed  in  the  Bousinesq  form,  that  is,  employing  eddy 
viscosity,  eddy  diffusion,  and  eddy  conductivity,  all  the  terms  in  the  equa¬ 
tions  are  time-averaged  quantities  and  no  need  exists  for  using  a  superscript 
bar  to  distinguish  the  time-averaged  quantities.  The  reader  is  cautioned  to 
remember,  however,  that  pv  represents  pv,  not  pv. 

The  equations  derived  are  essentially  identical  to  those  derived  in 

References  38  and  39  for  flow  over  a  flat  plate  with  uniform  free-stream 

0 

conditions,  except  that  the  diffusion  flux  used  here  is  in  a  more  general 
form  including  the  effects  of  thermal  diffusion  (Soret  effect)  and  unequal 
binary  diffusion  coefficients.  The  diffusion  introduced  by  pressure  gra¬ 
dients  and  body  forces  are  the  only  physical  phenomena  that  are  neglected. 

In  the  order  of  magnitude  arguments,  terms  of  the  following  type  have  been 
eliminated:  (1)  triple  cor’  -ilaticns  ,  (2)  derivatives  of  turbulent  correla¬ 

tions  parallel  to  the  wall,  and  (3)  correlations  involving  turbulent  com¬ 
ponents  of  molecular  transport  mechanisms. 
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2. 1.1.1  Conservation  of  specie  and  mass  equations 

A  mass  balance  of  an  individual  specie  in  a  unit  volume  results  in  the 
re lat ionship 


dx 


(pu.^) 


(pv^) 


dy  dy  ^  i^  +  f  l 


(2-1) 


where  y .  represents  the  rate  of  mass  generation  of  sp  ecies  i  per  unit 
volume  due  to  chemical  reactions,  in.  is  defined  in  terms  of  the  correlation 

ui 

of  the  fluctuating  components  of  concentration  and  normal  velocity,  that  is, 
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and  j;  is  the  mass  diffusion  rate  of  species  i  due  to  molecular  processes. 

This  latter  term  can  be  related  to  the  gradient  of  species  i  (as  well  as 

gradients  of  other  species),  the  family  of  binary  diffusion  coefficients, 

Df .  appropriate  to  the  boundary-layer  system,  and  the  multicomponent  thermal 
i  J  ij> 

diffusion  coefficients  through  the  Stef an-Maxwell  relations,  Equation 

8.1-3  of  Reference  39, 
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When  Equation  (2-1)  is  summed  over  all  the  species  in  the  system,  and 
it  is  remembered  that 


(2-4) 


because  of  the  definition  of  "m;  ss  diffusion"  and  that 


i 


because  of  a  conservation  of  atomic  species,  there  results 
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since  =  1.  Equation  (2-6)  is  the  familiar  overall  continuity  equation. 

When  Equation  (2-6)  is  considered  together  with  Equation  (2-1) ,  there 
results 
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2. 1.1. 2  Momentum  equation 

For  the  case  of  a  uniform  free-stream  condition,  the  momentum  equation 
can  be  written  as  follows: 
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where  the  eddy  viscosity  is  defined  in  terms  of  the  Reynolds  stresses  of 
turbulent  flow  by 


(Pv) 'u' 
(du/dy) 
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2. 1.1. 3  The  energy  equation 

The  energy  equation  for  this  general  system  can  be  written  as 
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where  the  DuFour  effect,  the  flux  of  heat  caused  by  thermal  diffusion,  is 
given  in  terms  of  the  diffusion  fluxes  of  the  pertinent  species  in  the  final 
term,  and  the  turbulent  enthalpy  transport  coefficient  is  defined  by 
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2.1.2  Binary,  non-reacting  system 

If  the  boundary  layer  is  constrained  to  a  binary  non-reu^uing  system 
(species  a  and  b) ,  considerable  simplification  results.  The  relations  for 
mass  transfer,  Equations  (2-7)  and  (2-3),  become 
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The  heat-transfer  equation  becomes 
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where  the  parameter  a  is  related  to  the  thermal  diffusion  coefficient  in 
Equations  8.1-17  and  8.1-10  of  Reference  39  by 
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and  is,  in  general,  relatively  invariant  with  concentration  (in  contrast  to 

T 

the  rather  large  variations  of  )  . 
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If,  additionally,  the  boundary  layer  is  constrained  to  the  cas^  of  con¬ 
stant  properties,  negligible  viscous  dissipation,  no  thermal  diffusion,  aid 
only  trace  amounts  of  foreign  gases,  the  final  simplified  two-dimensional 
relations  are  obtained.  These  are  for  mass  and  energy,  respectively 
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2.1.3  The  one-dimensional  mass,  momentum,  and  energy  equations 

The  general  two-dimensional  boundary-layer  relations  described  in  the 
preceding  subsections  are  such  a  complex  set  of  simultaneous,  non-linear, 
partial  differential  equations  that  their  complete  solution  including  even 
the  simplest  phenomenological  development  of  the  turbulent  transport  proper¬ 
ties  (tM,  tD,  and  e^)  has  yet  to  be  performed.  Many  of  the  complicating 
factors  in  this  equation  are  significant  only  in  the  region  close  to  the  wall. 

It  is  in  this  region  that  the  molecular  transport  properties  are  of  conse¬ 
quence  and  it  is  here  that  the  turbulent  transport  properties  are  most  variable. 
To  seek  a  two-dimensional  solution  in  this  region  seems  unwarranted,  particu¬ 
larly  in  light  of  the  extensive  empirical  evidence  that  it  is  characterized 
by  a  "law  of  the  wall."  Since  the  mathematical  statement  of  the  "law  of  the 
wall," 
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is  identical  with  the  implied  solution  of  a  one-dimensional  initial  value 
problem,  it  is  apparent  that  the  two  assumptions  are  compatible.  Thus,  it  is 
pertinent  to  look  at  the  relations  that  will  result  when  a  one-dimensional 
assumption  is  introduced. 

Tne  general  two-dimensional  continuity  equation  which  is 
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which  is  solved  to  yield 

,)V  =  (jv)  *  p  V 
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since  turbulent  correlations  vanish  at  the  surface,  (p'v1)  *  0.  Here,  p  v 

r  o  ’  ro  o 

can  be  considered  as  a  parameter  that  may  vary  with  x. 

The  one-dimens ional  assumption  leads  to  the  following  forms  for  the  mass, 
momentum,  and  energy  equations: 
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where  the  wall  shear  \  is.  in  general,  a  function  of  x,  as  is  q  which 
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where  ^injectant  t^le  ent^a^-Py  °f  the  injectant  evaluated  at  some  plane 

at  or  below  the  surface  and  ^SUppiied  crossing  the  same  plane  as  a 

consequence  of  conduction  or  radiation.  It  should  be  emphasized  that  hQ  is 
the  enthalpy  of  the  gas  mixture  adjacent  to  the  surface  as  contrasted  to  the 
enthalpy  of  the  pure  injectant. 

If  th~  boundary  layer  is  restricted  to  a  binary  nonreactive  system, 
these  equations  become 
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and  if  further  the  boundary  layer  is  constrained  to  an  essentially  one- 
component  system  containing  only  trace  amounts  of  foreign  species  and  having 
essentially  constant  properties,  the  final  relations  are  obtained. 
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The  last  three  relations  of  the  equations  above  can  be  expressed  in  terms 
of  a  new  set  of  variaoles  as  follows: 
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whore  the  viscous  dissipation  term  has  been  ignored  in  the  energy  relation. 
The  non-dimensional  variables  that  are  utilized  above  are  defined  by 
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In  addition, 

by 


the  momentum,  mass,  and  energy  transfer  coefficients  are  defined 
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The  solution  of  the  basic  equations,  either  one-  or  two-dimensional,  can 
be  divided  into  two  ma]or  aspects.  The  first  is  the  adequate  specification 
of  the  physical  parameters,  of  particular  importance  being  the  various  eddy 
transport  properties.  The  second  and  equally  important  aspect  involves  the 
mathematical  p  ocedures  for  solution  for  various  imposed  boundary  conditions. 

2.1.4  Int« gral  boundary -layer  equations  for  constant  fluid  properties 

Because  of  the  piysical  requirement  for  a  knowledge  of  the  turbulent 
transport  coefficients,  it  is  desirable  to  present  the  boundary-layer  equa¬ 
tions  in  a  fashion  permitting  the  experimental  evaluation  of  these  parameters. 
A  way  of  obtaining  these  parameters  from  mean  value  data  is  through  use  of 
integral  boundary-layer  equations.  These  equations  have  been  derived  for  the 
constant  property  boundary  layers  and  are  summarized  below.  These  relations 
are  presented  in  terms  of  the  internal  or  partial  values  (y  x.  5)  of  the  dis- 


placement,  momentum,  mass, 

and  enthalpy  thicknesses. 

These  thicknesses  are 

defined,  respectively,  as 

y 

r 

o*  -  /  (1  -  u)  dy 

o 

(2-44) 

y 

r 

y  =  /  u ( 1  -  u )  uy 

o 

(2-45) 

2-10 


v 

r  ud 


-  x )  ay 


(2-46) 


y 

r  u  ( i 


-  T)  dy 


(2-47) 


The  integral  relations  become  for  the  momentum  equation 
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for  the  energy  equation 
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and  for  the  mass  equation 
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When  y  ^  b  these  relations  reduce  to  the  more  conventional  integral  boundary- 
layer  equations.  These  are  as  follows  for  momentum,  energy,  and  mass  conser¬ 
vation,  respectively: 
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where  the  shape  factor,  H,  is  equal  to  the  ratio  of  the  displacement  thickness 
to  the  momentum  thickness  (&*/&). 

2.2  Solution  of  the  One -Dime  isional  Boundary-Layer  Equations 

In  this  section,  the  solutions  of  one-dimensional  boundary-layer  equations 
presented  in  the  preceding  section  will  be  developed.  These  solutions  are 
primarily  dependent  on  the  development  of  an  appropriate  model  of  turbulence 
and  its  coupling  with  the  one-dimensional  equations.  Certain  basic  considera¬ 
tions  will  be  discussed  with  regard  to  the  formulation  of  a  universal  "law  of 
the  wall."  Subsequently,  the  more  empirical  models  of  turbulence  will  be  con¬ 
sidered,  specifically,  the  pragmatic  mixing  length  concept.  A  functional 
dependence  of  the  mixing  length  on  the  various  boundary- layer  variables  will 
be  postulated  based  or.  certain  limiting  criteria  and  on  dimensional  arguments. 

This  postulated  mixing  length  relation  will  be  used  to  predict  profiles 
of  velocity,  temperature,  and  concentration  throughout  the  one-dimens ional 
portion  of  the  boundary  layer.  Extensive  re-evaluation  of  the  existing  data 
was  required  in  order  to  establish  the  significance  of  these  comparisons. 
Because  of  the  importance  of  these  evaluations,  they  will  be  described  in 
reasonable  detail. 

2.2.1  The  law  of  the  wall 

Mathematically  it  is  required  that  y  w  a  for  the  one-d.. mens  ional 
boundary-layer  equations  developed  in  the  preceding  section  to  be  strictly 
applicable.  In  practice,  however,  this  "one-dimensional"  regime  includes 
three  regions  of  interest  which  extend  over  a  significant  portion  of  the 
boundary-layer  thickness  as  indicated  in  Figure  1.  These  regions  include  the 
laminar  "sublayer"  wherein  molecular  transport  mechanisms  prevail,  the  tran¬ 
sitional  regime,  and  the  turbulent  core  of  the  velocity  profile  (the  log-law 
regime).  Inspection  of  the  basic  one-dimensional  equations  for  mass,  heat, 
or  momentum  transfer  indicates  that  a  solution  >_an  be  obtained  for  a  speci¬ 
fied  set  of  surface  conditions,  provided  only  that  the  variation  of  the  var¬ 
ious  thermodynamic  and  transport  properties,  both  molecular  and  turbulent, 
be  known.  The  turbulent  viscosity,  thermal  conductivity,  and  diffusion  coef¬ 
ficients  are  by  far  the  most  elusive  of  these  properties.  With  the  possible 
exception  of  these  turbulent  properties,  the  required  properties  are  functions 
of  local  conditions  (including  normal  fluxes).  Thus,  the  assumptions  implicit 
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in  the  acceptance  of  the  one-dimensional  equations  are  consistent  with  the 
empirically  accepted  "law  of  the  wall."  The  implication  of  the  "law"  is  that 
all  properties  of  the  boundary  layer  within  its  regime  are  functions  of  the 
normal  coordinate  from  the  wall  and  properties  at  the  wail,  for  example, 

T  -  ,;(y,  Tq,  qQ,  tq,  v  ,  physical  properties  of  mainstream  and  injected  gases) 

For  this  generalization  to  apply,  however,  it  is  necessary  that  the  turbulent 
transport  properties  be  the  function  of  the  same  set  of  variables.  It  should 
be  emphasized  that  there  is  no  obvious  analytical  logic  which  would  indicate 
this  simplified  dependence.  In  the  region  y  0,  the  equations  governing 
the  turbulent  components  of  velocity  cannot  be  reduced  to  a  one-dimensional 
set,  rather  a  three-dimensional  set  persists.  The  solution  of  this  set  requires 
initial  conditions  both  at  the  wall  and  at  two  mutually  perpendicular  planes 
normal  to  the  wall.  Considerable  data  exist,  however,  which  indicate  that  a 
"law  of  the  wall"  represents  the  mean  flow  and  that  properties  within  the  realm 
of  this  law  are  functions  only  of  the  normal  coordinate  and  the  wall  state. 
Clauser,  Reference  40,  was  one  of  the  first  to  make  extensive  use  of  this 
fact  under  conditions  of  severe  adverse  pressure  gradients .  The  implications 
of  these  facts  are  significant  in  that  they  indicate  the  direction  possible 
solutions  should  take. 

2. 2. 1.1  The  turbulent  transport  properties 

Because  of  the  current  lack  of  understanding  of  turbulent  mechanisms, 
"theoretical"  predictions  of  the  variation  of  turbulence  near  the  wall  must 
rely  on  empirical  input  into  relations  based  on  some  phenomenological  depend¬ 
ence.  Because  of  the  generality  of  the  ultimate  ^oals  of  this  analysis  and 
of  the  strong  desire  to  approximate  the  physical  situation,  certain  prerequi¬ 
sites  were  established  for  the  turbulent  transport  relations.  These  were: 

(a)  The  relations  must  indicate  a  continuous  variation  of  the  turbulent 
transport  properties  from  the  wall  to  the  fully  tuxbulent  region. 

(b)  The  relations  must  be  generally  applicable  to  mass,  momentum,  and 
energy  transport. 

(c)  The  relations  must  be  extendable  to  compressible  real  gas  property 

flow. 

( d )  The  relations  should  be  suitable  for  transpired  and  untranspired 
boundary  layers  without  any,  or  a  minimum,  modification  of  form. 

Two  basic  types  of  viscous  layer  hypotheses  have  been  proposed  in  the 
past.  The  first  type  predicts  the  variation  of  turbulent  viscosity  from  the 
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wall  to  the  fully  turbulent  regime.  Reichardt,  Rannie,  and  Deissler 
have,  in  References  41  through  43,  for  example,  proposed  such  varia¬ 
tions.  The  second  type  of  hypothesis  relates  to  a  variation  of  mixing  length 
from  the  wall  into  the  fully  turbulent  portion  of  the  boundary  layer.  Rotta, 
von  Karman,  and  van  Driest  (Refs.  44  to  46)  have  adopted  this  procedure.3 
Since  data  indicate  that  surface  mass  addition  strongly  affects  the  eddy 
viscosity  profile,  it  was  found  that  the  first  type  of  hypotheses  could  not. 
be  simply  modified  to  predict  this  variation.  On  the  other  hand,  the  success 
of  the  pragmatic  mixing  length  theory  in  predicting  profiles  in  the  fully 
turbulent  portion  of  the  boundary  layer  with  surface  mass  addition  has  been 
noted  (c.g.,  Refs.  13  and  17).  It  has  generally  been  concluded  that  the 
slope  of  the  linear  relation  between  mixing  length  and  distance  from  the  wall 
is  insensitive  to  surface  mass  addition.  As  a  consequence  of  this  apparent 
generality  of  the  mixing  length  approach,  it  was  adopted  for  the  present 
studies . 

2. 2. 1.2  The  mixing  length  postulate 

The  basic  mixing  length  postulate  can  be  expressed  as 


-  <-**  (§y  -  -v,  # 


(2-54) 


wher«  the  mixing  length,  I,  is  a  combination  of  various  correlations,  but 
retains  some  relationship  to  the  scale  of  turbulence.  Prandtl  proposed  that 
this  length  will,  in  its  simplest  form,  be  related  to  the  distance  from  a 
wall,  at  least  in  the  region  of  development  of  turbulence.  His  proposition 
that 
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3It  is,  of  course,  obvious  that  the  turbulent  viscosity  and  mixing  length 
postulates  are  simply  related  and  that  it  would  be  possible  to  cast  the 
former  in  terms  of  the  latter.  When  properties  are  constant  and  the  surface 
is  impervious,  the  relation  becomes 


i+  -  C—  (l  +  — 

v  \  €m 


or  when  £  >>> 

M 


vv 


l* ..  . 


2-14 


has  been  tested  under  a  variety  of  conditions  and  found  to  be  quite  adequate 
in  the  fully  turbulent  portion  of  the  law  of  the  wall  regime. 

As  the  wall  is  approached,  however,  this  simple  relation  is  no  longer 
appropriate  and,  in  fact,  it  can  be  shown  theoretically  that 
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This  is  a  consequence  of  the  Reichardt-Elrod  criterion  (see  Refs.  41  and  47, 
Appendix  A)  which  states  that 
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and  in  addition  that 
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if  the  flow  is  independent  of  the  spanwise  coordinate  and  if  axial  variations 
of  mean  flow  parameters  can  be  ignored  at  the  wall.  The  expression  for  the 
mixing  length  decay  (Eq.  (2-57))  is  required  by  Equation  (2-58)  but  also 
encompasses  Equation  (2-59)  . 

Thus,  two  criteria  have  been  specified,  namely,  Prandtl's  hypothesis 
which  will  be  taken  as  appropriate  in  the  fully  turbulent  portion  of  the 
"law  of  the  regime"  and  the  Reichardt-Elrod  criterion  as  expressed  by 
Equation  (2-57) . 

When  the  mixing  length  postulate,  Equation  (2-54),  is  introduced  into 
the  constant  property  one-dimensional  momentum  (Eq.  (2-29)),  there  results 
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where  f  +  -  iu^/v.  Equation  (2-60)  can  be  integrated  by  finite  difference 
techniques,  given  the  form  of  t+ ,  In  addition,  Equations  (2-28)  and  (2-30) 
can  be  simultaneously  integrated  with  Equation  (2-60)  to  yield  the  desired 
non-dimensional  concentration  and  temperature  distributions,  given  the  form 
of  l  ,  and  the  variation  (if  any)  of  turbulent  Prandtl  and  Schmidt  numbers. 

2.2.2  The  momentum  boundary  layer  without  surface  mass  addition 

Since  experimental  data  must  serve  as  the  principal  basis  of  the  appro¬ 
priate  mixing  length  relation  and  since  the  majority  of  the  available  boundary- 
layer  data  are  for  the  case  of  no  surface  mass  addition,  it  is  appropriate 
to  consider  this  simpler  case  first.  In  this  section  the  various  relations 
which  have  been  postulated  for  the  mixing  length  will  be  evaluated  in  light 
of  the  data  and  the  various  criteria  developed  in  the  preceding  section. 

2. 2. 2.1  Previous  mixing  length  representations 

Based  on  experimental  data,  several  functional  representations  of  the 
mixing  length  have  been  generated  which  satisfy  the  criterion  proposed  in 
Section  2.2.1,  along  with  others  that  do  not.  Some  of  the  better  known  rela¬ 
tions  are  included  in  Table  I.  This  table  also  indicates  the  ranges  of 
validity  for  the  various  expressions  and  their  compliance  (or  noncompliance) 
with  various  criteria. 

For  the  purposes  of  comparison,  the  emr  instants  for  each  of 

these  expressions  have  been  re-evaluated  consistent  with  the  following 
requirements : 

(a)  The  mixing  length  constant  k  0.44. 

(b)  A  value  of  u+  of  16.2  at  y+  ■  100. 

These  values  are  consistent  with  a  universal  log-law  velocity  profile  of  the 
form 

u+  -  5.2  ^ogiQy+  +  5.8  (2-61) 

where  the  constant  5.2  (fogg  10/fc ,  k  m  0.44)  has  been  selected  based  on 
evaluation  of  impervious  wall  velocity  profiles  from  a  variety  of  investiga¬ 
tors,  of  particular  significance  in  this  regard  being  the  data  of  Reference 
50  (see  Fig.  16  of  that  reference).  The  lowjr  values  of  a.  generally 
employed  have  been  derived  based  on  massed  data  plots  such  as  Figure  4  of 
Reference  40.  With  such  plots,  the  departures  of  u+  from  linearity  above 
the  curve  at  large  values  of  fog  y+  (the  wake  effect)  and  below  the  curve 
at  low  values  of  fog  y+  (Laminar  flow  effects)  tend  to  prejudice  the 
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selection  erf  the  linear  representations.  On  the  other  hand,  individual  pro¬ 
file  plots  (e.g.,  Fig.  5  of  Ret.  40)  tend  to  support  the  higher  value  of  \ 
used  here.  At  a  value  of  y  of  100  there  is  rather  uniform  agreement  amongst 
various  investigators  as  to  the  value  of  u+.  It  was  for  this  reason  that 
this  intercept  was  chosen  as  a  basis,  the  value  of  u+  of  16.2  being  a  repre¬ 
sentative  average  of  the  amassed  profile  data  (see  Refs.  3  and  40).  In 
Table  I,  six  flow  postulates  are  considered.  They  have  been  selected,  not 
with  an  eye  toward  completeness,  but  rather  in  an  effort  to  present  the  range 
of  postulates  available.  Historically,  the  two-layer  Prandtl-Tay lor  model 
must  be  token  as  the  point  of  departure  for  all  subsequent  models.  It  ta.es 
on  greater  significance  in  this  study  since  it  has  been  used  almost  exclu¬ 
sively  as  a  basis  for  past  mixing  length  studies  of  the  transpired  turbulent 
boundary  layer.  The  van  KfirmAn  three-layer  profile  prediction  was  introduced 
in  order  to  improve  the  u  relation  in  the  range  of  y  from  acout  5  to  30. 
In  so  doing,  however,  it  provides  no  particular  insight  into  the  mechanisms  in 
this  transitional  range.  Rotta,  by  perturbing  the  original  mixing  length 
expansion  to  account  for  a  constraining  effect  of  the  laminar  sublayer  and  by 
including  the  contribution  of  the  molecular  transport  properties  beyond  the 
sublayer,  achieved  continuity  of  shear  and  velocity  without  introducing  undue 
complication.  Reichardt  postulated  a  direct  expression  for  the  evaluation  of 
the  turbulent  viscosity  which  was  consistent  with  the  Reichardt-Elrod  cri¬ 
terion.  Similarly,  van  Driest  proposed  a  mixture  length  expression,  also 
consistent  with  this  criterion.  The  mixing  length  expression  utilized  in 
the  current  work,  and  developed  in  a  more  general  form  in  the  next  subsection, 
is  shov n  in  Table  I  in  its  reduced  form  for  the  case  of  constant  properties 
and  no  transpiration.  Additional  models  and  their  comparisons  are  discussed 
in  References  1  and  2.  For  each  model  given  in  Table  I,  the  requisite  mixing 
length  expression  is  presented,  if  the  original  representation  of  the  model 
was  in  some  alternate  form. 

As  would  be  expected,  the  three  multiple  layer  models  for  the  law  of  the 
wall  fail  to  satisfy  all  the  continuity  tests.  From  the  limiting  expression 
for  mixing  length  as  the  wall  is  approached,  it  can  be  argued  that  all  of  the 
models  given  satisfy  the  Reichardt-Elrod  criterion,  if  the  criterion  is 
limited  to  the  requirement  that  and  di+/dy+  approach  zero  at  the  wall, 

but  only  the  continuous  representations  do  so  in  a  non-degenerate  fashion. 

The  van  Driest  representation  of  f  +  includes  an  inflection  point  which  dis¬ 
tinguishes  it  from  the  other  models.  The  model,  to  be  developed  in  Section 
2. 2. 2. 2,  appears  to  combine  the  displacement  effect  of  Rotta  with  the  exponen¬ 
tial  decay  of  van  Driest. 
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This  latter  model  has  been  used  as  a  standard  of  comparison  in  the  formu¬ 
lation  of  Figure  3  wherein  the  differences  between  the  prediction  of  u+  by 
the  various  techniques  are  amplified.  It  is  apparent  that  the  three  single 
layer  models  agree  rather  well  and  any  selection  betveeen  them  should  be  based 
on  grounds  other  than  this  comparison. 


2. 2. 2. 2  Proposed  mixing  length  representation 

The  discussion  of  Section  2.2.1  leads  to  two  generalizations,  one  based 
on  experimental  evidence,  the  other  on  theoretical  considerations.  These  are 
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Several  means  of  expressing  a  relation  covering  the  full  range  of  y  and 
including  these  limiting  criteria  have  been  presented  in  Table  I.  Since 
physically  these  relations  are  concerned  with  the  transition  from  laminar  to 
turbulent  flow,  it  is  advantageous  in  considering  extensions  of  mixing  length 
theory  to  establish  some  physical  logic  for  the  selected  relation.  Unfortun¬ 
ately,  the  understanding  of  the  transitional  process  has  not  reached  a  state 
permitting  any  quantitive  specification.  Therefore,  the  selected  model  may 
be  based  only  on  qualitative  understanding  of  the  process,  dimensional  con¬ 
siderations,  and  the  above  limiting  criteria.  These  criteria  are  immediately 
satisfied  by  a  simple  implicit  relation  of  the  form 

3“  “  ('ey  ”  f)  (2-64) 


which  implies  that  the  rate  of  increase  of  the  mixture  length  is  proportional 
to  the  difference  between  value  postulated  by  Prandtl  (/cy)  and  its  actual 
value.  This  rate  of  increase  should  also  be  augmented  by  the  local  applied 
shear  and  retarded  by  the  local  viscosity.  Using  these  parameters  to  non- 
dimensionalize  the  above  relation  yields 

d7  *  •  £)~"J  (2-65) 

For  constant  shear  through  the  boundary  layer  (no  surface  mass  addition) , 
this  relation  becomes 


2-18 


where  y_  +  is  a  dimensionless  constant 

C4 

relation  yields 


which  for  large  y+  reduces  to  Rotta’s 


of  proportionality.  Integrating  this 


expression,  namely, 


(y*  -  O 


(2-68) 


Table  I  presents  Equation  (2-67)  along  with  five  of  the  more  conventional 
relations  previously  suggested.  In  addition,  this  expression  served  as  the 
basis  of  comparison  in  Figure  3  where  it  was  found  to  be  in  essential  agree¬ 
ment  with  other  continuous  models.  Because  of  its  conceptual  simplicity  and 
its  apparent  combination  of  the  better  features  of  Rotta's  and  van  Driest' s 
models,  Equation  (2-66)  has  been  chosen  as  the  basis  for  the  data  comparisons 
of  the  next  subsection.  In  general  form,  Equation  (2-65)  permits  inclusion  of 
the  change  in  shear  through  the  boundary  layer  introduced  by  surface  mass 
injection  in  a  physically  reasonable  manner. 


2. 2. 2. 3  Comparison  of  proposed  model  with  data 

The  adequacy  of  any  turbulent  boundary-layer  technique  is  measured  by 
the  degree  to  which  the  method  reproduces  the  experimental  data,  when  an 
accounting  is  made  of  the  uncertainties  which  may  exist  in  these  data.  The 
uncertainties  in  experimental  data  which  would  yield  errors  in  a  universal- 
velocity  profile  u  +  (y  +  )  are 

(a)  Uncertainties  in  friction  factor. 

(b)  Uncertainties  in  the  measured  local  velocities. 

It  is  common  knowledge  that  experimental  velocities  determined  by  pres¬ 
sure  probes  are  subject  to  errors  of  unknown  magnitude  when  the  probes  are 
in  the  near  vicinity  of  a  wall,  as  in  the  viscous  sublayer.  In  addition,  the 
effect  of  local  turbulence  can  be  significant  (as  will  be  noted  1  iter) .  To 
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some  extent,  errors  in  measured  velocities  using  hot-wire  anemometers  are 
correct ib 1 e ,  but  are  still  subject  to  some  uncertainty.  Friction  factors 
determined  from  the  inteyiaL  momentum  equation  (Eq.  (2-51)},  are  subject  to 
errors  primarily  because  they  require  the  presumption  of  two-dimensional  flow 
(no  spanwise  variations)  and  the  differentiation  of  experimental  data,  which 
amplifies  inherei  _  errors.  This  problem  is  greatly  magnified  for  transpira¬ 
tion  data,  as  will  be  seen  later,  in  that  the  friction  factor  is  interpreted 
by  a  small  difference  of  two  relatively  large  uncertain  numbers,  that  is, 
dO/dx  -  Vq/u^  For  the  purpose  of  comparing  the  various  mixing  length  tech¬ 
niques  to  experimental  data,  the  authors  have  chosen  the  following  data 
sources  as  subject  to  a  minimum  of  criticism*. 

(1)  Data  of  Laufer  (Ref.  48) 

These  data  were  achieved  for  fully  developed  turbulent  flow  in  pipes. 
Velocities  were  measured  by  hot-wire  anemometers,  and  friction  factors  were 
interpreted  from  the  measured  pressure  gradient.  The  major  uncertainty  here 
is  the  point  in  the  universal  velocity  distribution  where  flow  o.n  a  pipe 
departs  from  the  characteristics  of  a  flat  plate  and  the  origin  and  descrip¬ 
tion  of  the  wake. 

(2)  Data  of  Klebanoff  (Ref.  49) 

These  data  were  achieved  for  an  artificially  thickened  turbulent  boundary 

layer  on  a  flat  plate.  Velocities  were  measured  by  hot-wire  anemometers,  and 

turbulent  shears  were  measured  in  the  boundary  layer  near  the  wall  from  which 

the  surface  friction  factor  could  be  indirectly  determined,  that  is,  when 

v  -0,  3p/bx  *  0 
o 


for  small  y,  or 

t  *  t  ■  pu 1 v  ' 
w  r 

for  small  y  but  where  £  >>  v. 

m 

(3)  Data  of  Smith  and  Walker  (Ref.  50) 

These  data  were  measured  for  turbulent  boundary  layers  on  a  flat  plate. 
Velocities  were  determined  by  pressure  probes.  Friction  factoss  were  deter¬ 
mined  with  a  floating  element  flush  with  the  plate  surface. 


The  analytical,  universal  velocity  distribution  without  transpiration 
computed  using  Equations  (2-60)  and  (2-67)  is  compared  with  the  above  data 
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in  Fiaures  4(a),  4(b),  and  4(c).4  It  may  be  concluded  that  this  technique 
for  computing  velocity  distributions  does  a  good  job  of  representing  the 
available  data,  particularly  in  vie1*  of  the  data  of  Klebanoff,  and  those  of 
Smith  and  Walker. 

In  addition,  it  is  possible  to  compare  the  predictions  of  various  tur¬ 
bulence  parameters  involving  derivatives  of  the  velocity  profile  with  data. 
These  comparisons  are  presented  in  Figures  5  tnrough  7.  Fi  ire  5  compares 
the  predicted  turbulent  shear  variation  with  data  as  obtained  by  Laufer.  In 
Figure  6,  the  eddy  viscosity,  as  obtained  from  Schubauer 1 s  presentation 
vFig.  6  of  Ref.  51)  of  the  interpreted  data  of  Laufer  anu  Klebanoff,  are  com¬ 
pared  wi^h  the  predictions  of  the  model.  Finally,  in  Figure  7,  predictions 
of  two  r^.ated  functions,  the  energy  dissipation  function  and  the  turbulent 
energy  production  are  compared  with  the  values  interpreted  from  the  data  oi 
Laufer  and  Klebanoff  and  presented  by  Schubauer.  In  all  cases,  the  compari¬ 
sons  are  excellent. 

It  should  be  noted  that,  to  this  point,  comparisons  with  data  have  been 
based  primarily  on  surface  shear  data  obtained  in  some  direct  fashion  o  each 
station,  that  is,  turbulent  shear  measurements  by  floating  elements,  turbu¬ 
lence  corre laticnu ,  or  pressure  gradients  in  pipe  flow.  Thus,  the  uncertain¬ 
ties  recounted  at  the  beginning  cf  this  subsection  have  been  avoided  in  the 
verification  of  the  proposed  law— cf- the -wall  representation  for  the  case  where 
vq  -  0.  In  general,  this  verification  would  apply  equally  well  to  the  other 
continuous  models  compared  in  Figure  3,  Aeichardt  or  van  Driest. 

That  the  simple  model  proposed  in  Section  2. 2. 2. 2  accurately  predicts 
events  in  the  range  of  the  "law  of  the  wall"  seems  well-established  for  the 
constant  property,  impervious  wall  boundary  layer.  The  success  of  the  model 
when  extended  to  more  general  boundary  layers  will  be  considered  next. 

2.2.3  The  one-dimensional  momentum  boundary  layer  with  surface  mass 
addition 

Using  the  preceding  discussion  on  the  momentum  boundary  ayer  without 
surface  mass  addition  as  a  basis,  this  section  will  introduce  the  complicat¬ 
ing  factors  associated  with  surface  mass  addition.  In  the  past,  thi.->  exten¬ 
sion  has  been  performed  most  successfully  by  use  or  the  mixing  length  concept, 
its  success  having  led  to  its  adoption  in  the  present  study.  A  thorough  eval¬ 
uation  of  tne  existing  data  ruder  conditions  of  surface  mass  addition  will  be 

4Direct  comparison  of  the  flat-plate  relations  can  be  made  with  pipe  flow 
when  attention  is  confined  to  the  region  close  to  the  wall,  v/r  \  0.1. 
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presented  in  the  following  sections  along  with  the  general  comparison  of 
these  data  with  the  predictions  of  the  postulated  mixing  length  relations. 

2.  2. 3.1  Previous  mixing  length  representations 

The  fact  that  mixing-length  theory  is  readily  applicable  to  the  case  of 
transp  'ation,  as  previously  stated,  has  been  demonstrated,  for  example,  by 
Rubesin,  Kendall,  and  Mickley  and  Davis,  References  22,  17,  and  15.  Although 
these  authors  have  assumed  a  two-layer  approach  similar  to  the  Prandt 1-Taylor 
model  and  thus  do  not  satisfy  the  current  ground  rules  of  a  continuous  veloc¬ 
ity  gradient,  their  successes  require  some  consideration.  Therefore,  prior 
to  developing  a  model  consistent  with  the  specified  ground  rules  for  this 
study,  a  brief  review  of  the  previous  developments  will  be  presented  in  the 
following  paragraphs. 

For  the  so-called  laminar  sublayer  (e^v  *  0,  *  0) ,  the  solution  for 

the  momentum  equation  is 


+ 

u 


(2-b9) 


In  the  two-layer  model  this  solution  is  presumed  to  apply  over  the  range 

0  v  y  s  ya 


where  y^+  is  tht  non-dimensional  thickness  of  the  "laminar  sublayer"  and 
will,  in  general,  be  a  function  of  v  +,  the  mass-transfer  parameter. 

For  the  values  of  y  greater  than  ya  ,  the  quantity  €M/v  is  assumed 
to  be  much  greater  than  1,  allowing  Equation  (2-29)  to  be 
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whose  solution,  with  appropriate  limits  of  integration, 


yields 


(2  70) 


(2-71) 


2-22 


for 

+  + 

ya  K  y  < 


Several  postulates  have  been  made  wit  regard  to  the  variation  of  the 
laminar  sublayer  thickness,  y  +,  with  transp^^adon.  Rubesin  (Ref.  22) 

a 

assumed  two  types  of  conditions  from  which  two  forms  of  this  variation  could 
be  determined,  that  is,  postulate  1 

u  +y  +  =  constant  (2-72) 

d  a 


This  assumes  the  Reynolds  number  based  on  sublayer  thickness  and  the  inter¬ 
face  velocity  is  a  constant  and  is  analogous  to  a  transition  Reynolds  number 
This  assumption  yields,  from  Equation  (2-69) 


+ 


(2-73) 


where  the  subscript,  o,  on  y  +  refers  to  the  value  without  surface  mass 

d 

addition.  An  alternative  assumption,  postulate  2,  is 


which  yields 


u 


a 


+ 


constant 


(2-74) 


Van  DrifJt,  Reference  2,  assumed  that  the  Reynolus  number  based  on  the 
thickness  of  the  laminar  sublayer  and  the  local  shear  velocity  was  invariant 
that  is, 


2-23 


or 


y , 


+ 


o 


1 


(2-76) 


Kendall  (Ref.  17)  found  that  the  variation  of  the  laminar  sublayer  thick¬ 
ness  required  for  Equation  (2-71)  to  fit  his  data  in  the  turbulent  core  was 
remarkably  similar  to  that  postulated  by  van  Driest  (Eq.  (2-76))  and  the 
second  postulate  of  Rubesin,  Equation  (2-74) .  Both  of  these  theoretical 
postulates  yielded  essentially  identical  values  of  laminar  sublaver  thickness 
as  indicated  in  Figure  8. 


2. 2. 3. 2  Proposed  mixing  length  representation  and  resultant  predictions 

Since  the  above  relations  depend  on  the  existence  of  a  purely  laminar 
sublayer  up  to  y  +  to  form  the  relation  for  y  +,  their  extension  to  the 

d  d 

techniques  permitting  a  continuous  decay  of  turbulence  all  the  way  to  the 
wall  is  not  well  defined.  Fortunately,  the  exteiu  on  of  the  technique  pro¬ 
posed  in  Equation  (2-65)  is  quite  straightforward.  This  relation  becomes 


dy 
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.j/j/iL 

VVPo 


(2-77) 


where  the  effect  of  transpiration  manifests  itself  in  the  locally  variable 
shear  ratio.  With  Equation  (2-26),  this  expression  becomes,  for  the  constant 
property  boundary  layer 


S ill 
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(2-78) 


Because  of  the  implicit  nature  of  Equation  (2-78) ,  it  is  intractable  in 
closed  form  for  vq+  f  0.  Since  the  use  of  a  computer  had  long  since  been 
accepted,  this  did  not  result  in  any  significant  complication  of  the  solu¬ 
tion  procedure. 

The  results  obtained  when  Equation  (2-78)  is  combined  with  Equation  (2-60) 
and  the  indicated  integration  is  carried  forward  numerically  are  presented  in 
Figure  9  for  a  range  of  values  of  the  blowing  parameter  vq+  and  yQ+  -  11.83, 
its  value  for  the  zero  surface  mass-transfer  case.  Before  an  assessment  can 


be  made  of  the  accuracy  of  this  inner  region  prediction  technique,  including 
the  assumption  of  a  constant  y  +,  it  will  be  necessary  to  demonstrate  means 

a 

of  obtaining  the  local  skin  friction  to  be  used  in  the  dimensionless  param¬ 
eters  uf  and  y  .  This  is  the  topic  of  the  next  section. 

2. 2. 3. 3  Evaluation  o£  available  friction  factor  data 

As  was  noted  in  the  preceding  section,  values  of  the  local  wall  shear 
cure  required  to  normalize  velocity  profile  data  if  these  data  are  to  be  com¬ 
pared  directly  with  the  predictions  of  the  law  of  the  wall.  Most  of  the 
existing  skin-friction  data  for  cases  with  transpiration  have  been  evaluated 
employing  the  integral  momentum  equation  for  two-dimensional  flow,  namely: 

Cr  ,  V  Q  du 

+  (2  +H)  — ±.  (2-79) 

2  dx  ux  ui  dx 

wherein  the  sum  and  difference  between  experimentally  determined  parameters 
are  required.  This  is  rather  an  inaccurate  procedure  when  the  positive  and 
negative  terms  are  both  large  and  of  the  same  magnitude.  Although  the 
accuracy  of  friction  factors  interpreted  from  data  by  use  of  Equation  ^ 2  — 79 ) 
is  reduced  when  vQ/ui  is  positive,  it  is  enhanced  when  v Q/^l  is  negative, 
that  is,  for  the  cases  of  suction.  Thus,  the  suction  data  of  Dutton,  Refer¬ 
ence  16,  becomes  an  important  basis  of  comparison. 

For  positive  values  of  v Q/'u1,  the  data  obtained  at  MIT  under  the  direc¬ 
tion  of  H.  S.  Mickley  represent  an  extensive  source  of  two-dimensional  data. 
Before  comparisons  are  made  with  these  data,  however,  it  is  well  to  subject 
it  to  some  scrutiny,  as  is  done  in  the  next  paragraphs. 

A  series  of  reports  and  doctoral  dissertations  containing  transpiration 
data  have  been  originated  at  MIT,  authored  by  Mickley,  Stewart,  Squyres , 
Gessner,  Curl,  Davis,  Kendall,  Smith,  Goodwin,  and  Butensky  (Refs.  9  through 
15  and  17  through  20).  These  sources  contain  data  for  heat,  mass,  and  momen¬ 
tum  transfer  under  the  influence  of  transpiration  into  incompressible  turbu¬ 
lent  boundary  layers  on  a  flat  porous  plate.  With  one  exception,  all  reported 
friction  factors  were  determined  from  velocity  profiles  and  the  integral 
momentum  equation.  In  detail,  uncertainties  in  this  method  of  evaluation  of 
the  friction  factors  arise  for  the  following  reasons: 

(1)  Lack  of  knowledge  of  local  transpiration  rates  in  the  vicinity  of 
the  survey  probe,  in  contrast  to  the  reported  average  compartment  transpiration 
rates . 
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(2)  The  influence  of  three-dimensional  wind-tunnel  flow  effects  when 
evaluations  are  based  on  the  presumption  of  two-dimensional  flow. 

(3)  As  blowing  rate  increases,  the  gradient  in  the  momentum  thickness 
(dc/dx)  approaches  the  blowing-rate  parameter,  v  /u, ,  and  the  calculated 
friction  factor  is  the  result  of  a  small  difference  in  two  large  numbers, 
each  with  some  error,  Equation  (2-79). 

(4)  The  determination  of  the  gradient  of  the  momentum -thickness  along 
the  plate  requires  differentiation  of  experimental  data. 

For  these  reasons,  an  attempt  has  been  made  here  to  report  f  momentum 
data  in  an  alternate  form,  in  order  to  distinguish  differences  in  experimen¬ 
tal  results,  and  in  the  hope  of  determining  experimental  friction  factors  in 
light  of  all  of  the  available  data.  To  this  end,  the  following  format  for 
data  presentation  was  adopted.  The  momentum  integral,  Equation  (2-79),  when 
integrated  with  respect  to  x,  yields 


Piui  f 

11  J 


x 


b  (2  +  H) 


d  in  uA 
dx 


dx  (2-80) 


The  terms  on  the  right  side  of  the  above  equation  can  be  evaluated  directly 
from  velocity  profile  data,  th?  surface  injection  rate  distribution,  and  the 
velocity  gradient  along  the  surface.  No  differentiation  of  experimental 
profile  data  is  required  and  the  attendant  uncertainties  of  such  an  operation 
are  avoided  in  the  comparisons  that  follow.  Although  (Cf/2)Rex  cannot  be 
considered  a  universal  relationship  independent  of  such  factors  as  the  loca¬ 
tion  of  the  boundary- layer  transition,  the  strength  of  boundary-layer  trips, 
and  the  streamwise  pressure  gradients,  this  relationship  is  quite  useful  in 
the  comparison  of  data  all  taken  with  similar  equipment  and  flow  conditions. 

The  drag  parameter  (C^/2)Rex  for  some  of  the  MIT  data  is  presented 
in  Figures  10(a)  through  (d) .  It  is  seen  that  a  significant  spread  in  the 
data  exists  between  and  within  the  data  sources.  This  spread  is  likely  to 
be  due,  primarily,  to  the  uncertainty  in  the  local  blowing  rate,  that  is,  in 
the  blowing  rate  along  the  axial  plane  in  which  the  velocity  profiles  were 
measured,  in  contrast  to  the  spa. wise  average  blowing  rate  that  was  fne  actual 
quantity  measured.  For  cases  in  which  concentration  or  temperature  profiles 
were  measured  simultaneously  (or  at  least  for  the  srr.e  operating  conditions) 
with  the  velocity  profiles,  mass  or  energy  balances  can  be  performed  on  the 
system  in  order  to  establish  the  consistency  of  the  reported  blowing  rates. 


2-26 


For  cases  in  which  concentration  profiles  were  measured,  the  following 
overall  mass  balance  can  written  on  a  given  specie: 

x  o 

/  vo  (%  -  xi1)  Jx  ■  /  u  (xi 
o  o 

where  x^2  is  the  mole  fraction  of  the  species  i  in  the  undiluted  injectant 
and  Xj^  is  its  mole  fraction  in  the  free  stream.  Equation  (2-81)  presumes 
that  the  mole  density  of  this  mixture  is  essentially  constant  (i.e.,  tempera¬ 
ture  and  pressure  are  constant  for  a  mixture  of  perfect  gases )b.  In  terms 
of  the  mass  thickness  (Eq.  (2-46)  at  y  =  5)  this  becomes 
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since  x^  and  Xj.  were  constant  with  distance  along  the  plate  for  the 
available  data. 

Similarly,  an  overall  energy  balance  in  low  speed  flow  with  constant 

C  results  in 
P 

1  r  *  r 

-  V  -  ^  J  dx  -  ^  J  -  TJdx  {2-83 ) 

P  o 

where  2q  ■  sum  of  electrical  heat  input  less  losses  per  unit  surface  area, 
and  T2  is  the  temperature  of  the  pure  injectant  prior  to  its  contact  with 
the  surface  heating  units. 

The  right-hand  side  of  Equation  (2-82)  was  plotted  against  its  left- 
hand  side  (Ref.  17),  usinc  the  helium  data  of  that  reference.  The  resultant 
curve  has  been  duplicated  here  as  Figure  11.  The  slope  of  the  line  best 
fitting  these  data  is  1.10.  Thus,  if  all  other  parameters  of  Equation  (2-82) 
are  presumed  accurate,  it  is  possible  to  conclude  that 
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(2-82) 


5For  the  available  data  this  assumption  is  less  restrictive  than  the  constant 
mass  density  assumption  introduced  in  the  theoretical  development. 
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for  cases  when  vq  is  approximately  constant.  The  required  presumption  is 
relatively  sound  since  no  differentiation  or  extrapolation  is  required  to 
obtain  the  ne  '■'ssary  parameters.  For  this  reason  the  blowing  rates  reported 
in  Reference  17  include  this  correction  factor. 

In  Figure  12(a),  the  right-hand  side  of  Equation  (<i-t3)  was  plotted 
against  its  left-hand  side  for  one  run  of  Squyres  (Refs.  9  and  10).  When 
heat  losses  are  ignored  it  is  seen  that  the  slope  of  the  curve  is  1.38  indi¬ 
cating,  for  approximately  constant  v  , 

v 

v-°«*MflU.  .  1>38 

°measured 


Equation  (2-83)  indicates  that  inclusion  of  heat  losses  would  further  increase 
this  factor.  In  Figure  12(b),  a  similar  plot  was  prepared  for  one  run  of 
Reference  12.  It  is  seen  that  a  fairly  good  energy  balance  is  obtained  for 
these  latter  data.  The  effect  of  heat  losses  is  also  indicated  based  on  the 
original  estimates  of  Curl  in  Reference  12.  He  considers  the  loss  calculation 
to  be  accurate  only  to  within  +50  percent,  however.  The  indicated  values  of 
the  ratio  v0acfcual  to  v0n)eagured  are  0.89  and  0.95  for  the  balances, 
excluding  and  including  heat  losses,  respectively. 

It  might  appear  unusual  that  this  parameter  varies  so  significantly 
from  one  set  of  data  to  another,  since  nonuniformities  in  blowing  distribu¬ 
tion  and/or  three-dimensional  effects  would  normally  be  considered  a  char¬ 
acteristic  of  the  apparatus.  In  each  instance,  however,  changes  were  made 
between  individual  experiments  in  the  configuration  of  the  porous  wall  and 
the  adjacent  side  walls. 

As  a  consequence  of  these  determinations,  the  transpiration  rates  of 
Squyres  were  "corrected"  for  the  particular  run  studied.  The  consequent 
values  of  Cf  Re^  were  also  modified  and  replotted  in  Figure  13.  There 
results  a  significant  improvement  in  the  correspondence  of  these  data  with 
those  of  the  other  investigators.  This  lends  credance  to  the  postulated 
error  mechanism. 

It  is  apparent  that,  with  the  exception  of  the  no-blowing  data,  a 
straight  line  is  the  only  justifiable  curve  through  the  assessed  data.  The 
resultant  relation!:  will  be  of  the  form 
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(2-84) 
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yielding  on  differentiation  of  the  entire  left  side 

C  . 

•f  “  ab  RexD_1  (2-85) 

The  following  table  summarizes  the  results  obtained  using  the  straight  lines 
indicated  in  Figures  10(a),  (b) ,  (d) ,  (e),  and  13: 


Vv  i 

a 

b 

ab 

b-i 

0.000 

5 . 50x10"** 

0.768 

4.23x10"^ 

-0.232 

.001 

3.26X10"2 

.784 

2 . 55xl0"2 

-  .216 

.002 

8.24X10"2 

.689 

5.67X10"2 

-  .311 

.003 

2.21X10"2 

.762 

1.68X10"2 

-  .238 

.004 

2.26X10"1 

.595 

1.34X10"1 

-  .405 

Although  a  number  of  tests  have  been  performed  at  nominal  transpiration 
ratios  of  0.005,  the  data  when  presented  in  the  form  of  Figure  10  scatter 
wildly.  As  an  example  of  this  difficulty,  one  run  from  Reference  17  is 
plotted  in  Figure  14.  It  is  obvious  that  the  evaluation  of  a  slope  from 
these  data  is  impossible.  It  should  be  emphasized  that  the  magnitude  of  the 
ordinate  of  the  plotted  points  represent  on  the  average  only  about  3  percent 
of  the  expeiimentally  determined  Reg. 

When  velocity  gradients  exist,  another  integral  form  of  the  momentum 
equation  is  more  convenient  than  Equation  (2-81)  in  the  evaluation  of  data, 
namely : 


H)(9u7  1+H  dux 

(2-86) 

where  the  final  term  can  usually  be  ignored  if  the  shape  factor,  H,  does  not 
vary  much  from  an  average  value,  H.  This  relation  was  used  with  two  runs  of 
Davis,  specifically  those  runs  which  he  used  to  define  the  friction  factor 
at  a  Vui  of  ^.005.  The  results  are  shown  in  Figure  15.  Again  the  same 
scatter  is  noticed,  and  any  conclusion  with  regard  to  actual  values  of  fric¬ 
tion  factor  must  be  seriously  questioned. 

Tewfik  (Ref.  35)  also  used  the  momentum  equation  to  evaluate  friction 
factors  at  transpiration  ratios  to  0.003.  His  determinations  indicated  the 
same  difficulties  described  here  for  the  MIT  data.  It  seems  apparent  that 
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these  determinations  have  also  been  affected  by  feedback  of  downstream  per¬ 
turbations  of  the  flow  field.  Such  perturbations  have  been  found  significant 
in  other  experiments  (Ref.  14).  For  this  reason,  the  values  of  skin  friction 
determined  in  the  midrange  of  his  test  section  are  considered  most  reliable. 

It  is  thus  apparent  that  the  values  of  skin  friction,  determined  by  the 
momentum  equation,  for  transpiration  ratios,  v^/ui#  in  excess  of  0.003  are 
subject  to  considerable  uncertainty  utilizing  boundary-layer  integrals. 

Because  of  the  above  uncertainties,  an  alternate  method  of  skin-friction 
evaluation  was  attempted  in  Reference  17.  It  was  based  on  interpretations 
of  the  dynamic  pressure  measurements  cf  the  pitot  probe  in  the  vicinity  of 
the  wall.  Basically,  two  assumptions  were  made.  The  first  was  that 


dy 


undisturbed 


(2-87) 


where  5^  was  the  probe  height  and  y  was  the  distance  of  the  midDlane  of 
the  probe  from  the  wall.  The  second  presumption  was  that  the  velocity  pro¬ 
file  ultimately  approached  the  wall  according  to  the  laminar  wall  law  rela¬ 
tion  (Eq.  (2-69)).  The  technique  is  described  in  detail  in  Appendix  B. 

It  is  concluded  in  Appendix  B  that  the  procedure  does  indeed  yield  valid 
friction  factors  and  that  the  accuracy  of  the  technique  does  not  degenerate 
as  the  mass  addition  rate  is  increased.  It  should  be  emphasized  that  the 
procedure  is  independent  of  any  assumption  with  regard  to  the  growth  (or 
decay)  or  turbulence  near  the  wall  except  that  reasonably  monatonic  behavior 
is  anticipated.  On  the  contrary,  such  profile  analysis  procedures  as  those 
proposed  by  Clauser  (Ref.  40)  require  the  pre-knowledge  of  the  "law  of  the 
wall"  turbulent  velocity  distribution.  Since  the  law  of  the  wall  is  signifi¬ 
cantly  affected  by  surface  mass  addition,  the  application  of  the  Clauser  pro¬ 
cedure  here  will  be  valid  only  where  the  law  of  the  wall  with  surface  mass 
addition  is  verified  by  such  independent  friction-factor  determinations  as 
were  discussed  in  this  subsection. 

In  Figure  16  the  friction  factors  developed  in  Reference  17  by  use  of 
the  technique  described  in  Appendix  B  are  presented,  together  with  those 
values  obtained  from  the  momentum  equation  (see  the  table  following  Eq.  (2-85)). 
In  addition,  a  curve  based  on  the  data  of  Smith  and  Walker  (Ref.  50)  is  also 
shown  along  with  the  Schultz-Grunow  correlation  curve.  The  general  agreement 
is  well  within  that  expected.  In  order  to  present  the  results  of  these 
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evaluations  in  a  more  useful  fashion,  Figure  17  was  prepared  based  on  a 
sequence  of  cross  plots  of  the  data  of  Figure  16.  This  latter  figure  might 
be  considered  the  "conclusion"  of  tnis  evaluation  of  available  friction- 
factor  data  for  the  turbulent  boundary  layer  with  surface  mass  addition. 
Values  will  be  taken  from  this  figure  in  the  next  section  to  formulate  the 
nondimensional  velocity  profile  parameters  and  their  subsequent  comparison 
with  the  proposed  law  of  the  wall. 

2. 2. 3. 4  Comparison  of  proposed  model  with  data 

Probably  the  most  meaningful  test  of  the  proposed  law  of  the  wall  rela¬ 
tions  is  its  "direct"  comparison  with  measured  profile  data.  Tnis  comparison 
had  to  await  the  prior  establishment  of  the  friction  factors  appropriate  to 
these  measured  profiles.  This  was  done  in  the  preceding  subsection  and 
resulted  in  a  generalized  plot  of  friction  factor  as  a  function  of  the 
Reynolds  number  and  the  blowing  ratio  {Fig.  17).  In  the  comparisons  to 
follow,  friction  factors  will  be  obtained  from  tnis  figure  for  all  flat- 
plate  profiles. 

A  series  of  comparisons  are  presented  in  Figures  18(a)  and  (b)  using 
the  blowing  data  of  Kendall  and  Gessner  and  the  suction  data  of  Dutton.  For 
the  suction  profiles,  essentially  asymptotic  condition  prevailed  and  thus 


assuring  rather  good  accuracy  in  the  determination  of  friction  factors  for 
these  data.  The  quality  of  all  the  comparisons  is  seen  to  be  good.  The 
deviations  noted  in  the  range  of  y+  up  to  30  are  most  likely  due  to  probe 
and  turbulence  corrections,  the  former  effect  being  significant  only  for  y+ 
less  than  10  (see  Fig.  B-4) .  Turbulence  contributions  to  the  dynamic  pres¬ 
sure  measurements  cannot  be  evaluated  without  introducing  multiple  assump¬ 
tions  but  appear  to  approximate  the  noted  discrepancy  (based  on  the  u'u' 
correlation  data  of  Refs.  48,  49,  and  19). e 

6 A  mixing  length  model  was  proposed,  namely 


where 


u'u’ 


(dllf 

Vdy  J 


Using  this  model  good  agreement  was  obtained  with  all  available  data. 
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Certain  additional  data  have  been  reported  by  Tewfik  (Ref.  35) .  The 
woven  screen  used  in  these  experiments  produced  definite  roughness  efforts. 

The  electromesh  screen  used  at  MIT  did  not  produce  these  effects  and  thus 
direct  comparison  is  difficult.  Using  the  data  obtained  without  blowing  it 
was  noted  that  the  required  value  of  y  +  for  the  proposed  model  (Eq.  (2-of)) 

d 

was  6.86  (as  compared  to  11.825).  Using  this  value,  a  comparison  was  made 
with  one  of  the  profiles  reported  in  Reference  35,  the  results  being  shown 
in  Figure  19.  The  comparison  is  quite  satisfactory.  The  reduction  of  the 
effective  laminar  sublayer  thickness  by  wall  roughness  is  in  accord  with  the 
supposition  of  Rotta.  It  is  significant,  however,  that  blowing  does  not 
appear  to  lessen  this  roughness  effect. 

Final  conclusions  with  regard  to  the  adequacy  of  the  proposed  wall  law 
model  will  be  reserved  until  after  concentration  and  temperature  profile 
comparisons  are  presented.  To  this  point  a  certain  internal  consistency  has 
been  established  between  the  chosen  model,  the  measured  profiles,  and  the 
evaluated  friction-f  i.ctor  data.  The  weak  link  is  of  course  the  accuracy  of 
the  friction-factor  data.  Mass-transfer  and  heat-transfer  coefficients  are 
more  directly  determined  and  thus  will  serve  as  a  better  basis  for  the  final 
assessment  of  the  model. 

2.2.4  The  one-dimensional  diffusion  and  thermal  boundary-layer  equations 
with  surface  mass  addition 

2.2.4. 1  Mixing  length  form  of  equations 

When  viscous  dissipation  and  chemical  reactions  do  not  occur,  and  when 
the  foreign  gas  is  present  in  trace  quantities  so  that  properties  are  essen¬ 
tially  constant  across  the  boundary  layer,  the  diffusion  and  energy  equations 
in  the  inner  portion  of  the  boundary  layer  can  be  written  as  fellows: 


dy  + 


axl.  !ili JiiaXl 

dyt  1  + 1^- 1*‘  *4 


(2-88) 


(2-89) 
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In  these  equations,  the  mixing  length  term  and  the  velocity  gradient  are 
ootained  in  numerical  form  from  previous  numerical  integrations  of  Equations 
(2-78)  and  (2-60).  Integration  of  Equations  (2-88)  and  (2-89)  can  proceed 
numerically  once  values  of  the  turbulent  Prandtl  and  Schmidt  numbers  (Pr' 
and  Sc‘ )  are  specified. 

The  MIT  data  of  Curl  and  Kendall  were  used  to  evaluate  the  turbulent 
Prandtl  and  Schmidt  numbers  employed  in  the  numerical  work  of  this  study. 

If  in  the  above  diffusion  equation  it  is  assumed  that  the  turbulent  eddy 
transport  ie  much  larger  than  the  molecular  transport  and  vq+  is  set  equal 
to  zero,  it  can  be  shown  readily  that 

x+  -  (Sc')u+  +  const  (2-90) 

In  Kendall's  experiment,  Reference  17,  he  introduced  a  helxum-air  mixture 
into  an  air  boundary  layer  at  such  small  mass  injection  rates  that  an  equa¬ 
tion  of  the  form  of  Equation  (2-90)  would  be  expected  to  hold  over  a  portion 
of  the  boundary  layer.  When  hi.s  concentration  data  are  plotted  against  thv 
local  dimensionless  velocity  as  in  Figure  20,  it  is  found  that  a  straight- 
line  relationship  holds  quite  well  for  the  data  and  that  the  turbulent 
Schmidt  number  can  be  represented  by  0.75.  The  dimensionless  velocity  range 
used  here  is  such  that  molecular  viscosity  effects  can  be  neglected  and,  yet, 
the  one-dimensional  character  of  the  flow  should  still  apply.  Under  condi¬ 
tions  of  blowing,  where  one -dimensional  eddy  transport  prevails,  the  turbu¬ 
lent  Schmidt  number  can  be  expressed  as 


Sc'  - 


1  +  v  V  .  + 

_ S _ _  &L_ 

1  +  v  x  du 
o 


(2-91) 


i 


Kendall  found  that  he  could  not  evaluate  his  data  with  blowing  as  accurately 
as  the  very  small  (essentially  zero)  blowing  data  described  above;  therefore, 
in  the  present  work,  the  small  increase  in  turbulent  Schmidt  number  exhibited 
by  the  da^a  at  the  increased  blowing  rates  will  be  ignored.  Turbulent  Prandtl 
numbers,  in  the  neighborhood  of  the  wall,  were  found  from  the  temperatuie 
distribution  date  of  Curl  ar.d  Gessner  in  a  manner  similar  to  that  shown  in 
Figure  20.  It  was  found  that  a  Prandtl  number  of  0.75  fitted  the  data  well. 
Again,  under  the  conditions  of  transpiration,  these  data  contained  serious 
uncertainties  in  friction  factor  and  heat-transfer  coefficients  and  in  the 
velocity  and  temperature  distributions  in  the  inner  portions  of  the  boundary 
layer  where  these  one-dimensional  equations  can  be  expected  to  apply. 
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For  the  purposes  of  this  study,  therefore,  the  turbulent  Prandtl  and 
Schmidt  numbers  were  assumed  constant,  independent  of  surface  mass  transfer, 
and  each  equal  to  0.75.  These  simplifications  facilitated  the  integration 
of  the  diffusion  and  energy  equations,  the  results  of  which  are  shown  in  the 
next  section. 

2. 2. 4. 2  Comparison  of  concentration  and  temperature  distributions  with  data  1 

The  theoretical  parametric  predictions  of  the  law  of  the  wall  nonaimen- 
sional  concentration  (Sc  ■  0.21,  Sc1  «  0.75)  and  temperature  (Pr  ■  0.7, 

Pr'  *  0.75)  profiles  are  presented  in  Figures  21  and  22,  respectively.  The 
significance  of  these  figurec  will  become  apparent  in  the  data  comparisons 
that  follow  and  in  the  later  stages  of  the  analysis. 

The  theoretical  predictions  corresponding  to  the  experimental  concentra¬ 
tion  profiles  of  Kendall  (using  his  reported  friction  and  mass-transfer  fac¬ 
tors)  are  compared  with  his  data  in  Figure  23.  The  departure  of  the  experi¬ 
mental  data  from  the  theory,  in  its  range  of  applicability,  is  well  within  the 
uncertainty  of  Kendall's  reported  mass-transfer  factors;  that  is,  adjustment 
of  the  reported  mass -transfer  factors  of  less  than  10  percent  is  required  for 
the  predicted  distributions  to  agree  with  the  normalize!  experimental  concen¬ 
tration  profiles  for  the  three  blowing  rates  shown.  Note  that  the  wall  pres¬ 
sure  profile  technique,  discussed  in  Appendix  B,  was  used  to  determine  the 
experimental  friction  factors  needed  in  these  dimensionless  coordinates. 

The  theoretical  predictions  corresponding  to  the  temperature  distribution 
data  of  Curl  are  shown,  along  with  the  data,  in  Figure  24.  The  ordinate  in 
this  figure  is  the  dimensionless  temperature  parameter  defined  by  Equation 
(2-35) .  It  will  be  noted  that  two  sets  of  data  are  given  for  largest  blowing 

rate  where  v  +  «  0.168.  The  differences  in  these  data  are  the  Stanton  num- 
o 

bers  (Cjj)  employed  in  the  correlations;  the  data  indicated  by  square  symbols 
are  based  on  the  directly  measured  Stanton  numbers  as  reported,  and  the  data 
represented  by  diamond  symbols  are  based  un  Stanton  numbers  with  the  heat- 
loss  correction  excluded.  The  better  agreement  of  the  latter  representation 
should  not  be  considered  significant,  rather  the  difference  between  the  two 
curves  should  be  considered  as  one  measure  of  the  effect  of  the  uncertainties 
in  the  heat-transfer  coefficient  data.  It  should  be  noted  that  the  direct 
determination  of  heat-transfer  coefficient  is  linearly  dependent  on  the  local 
transpiration  ratio.  This  potentially  represents  the  most  significant  source 
of  uncertainty.  The  increasing  deviations  of  the  data  from  the  predictions 
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as  the  wall  is  approached  are  no  doubt  due  to  the  relatively  large  gradients 
of  temperature  close  to  the  wall  and  other  wall-probe  interactions. 

It  can  be  concluded  from  the  above  comparisons  of  data  and  theory  that 
th‘"  law-of-the-wall  predictions  of  temperature  distributions  (employing 
Pr  ■  0.7  and  Pr‘  ■  0.75)  and  helium  concentration  distribution  (using  i»c  *  0.21 
and  Sc'  -  0.75)  are  quite  satisfactory  when  consideration  is  given  to  the 
possible  uncertainties  in  the  data. 

2.2.5  Conclusions  regarding  law  of  the  wall 

A  law  of  the  wall  was  developed  in  the  previous  sections  that  produced 
continuous  profiles  of  velocity  in  the  laminar,  transitional,  and  fully  tur¬ 
bulent  regimes  near  the  wall  under  the  conditions  of  nonzero  surface  mass 
transfer  and  a  constant  deriitv  in  the  boundary  layer.  Profiles  of  the  con¬ 
centration  of  trace  elements  and  sma^l  temperature  variations  were  also 
evaluated.  The  theoretical  model  em  ..oyc  a  mixing-leng  \  formulation  that- 
continuously  approaches  the  Prandtl  requirement,  di/dy  -*  K,  in  the  fully  tur¬ 
bulent  region  and  the  Elrod  criteria  at  the  wall.  In  addition,  the  mixing- 
length  formulation  contain  factors  t«.at  can  be  treated  as  constants,  indepen¬ 
dent  of  the  surface  blowing  rate. 

When  the  theoretical  model  is  compared  to  velocity  profile  data,  after 
skin-friction  factors  have  been  established  to  permit  expressing  the  profiles 
in  dimensionless  form  characteristic  of  a  wall  law,  it  is  found  that  there  .s 
rather  good  agreement  over  a  considerable  range  of  surface  mass  transfer. 
Further,  temperature  and  concentration  profiles  were  predicted  and  compared 
with  data,  after  turbulent  Prandtl  and  Schmidt  were  evaluated  from  zero-blowing 
data  and  assumed  unchanged  by  surface  nass  transfer.  The  good  agreement  of 
the  latter  comparisons  further  established  the  usefulness  of  the  currently 
proposed  model,  at  least  in  the  vicinity  of  the  wall. 

As  stated  earlier,  wall  laws  are  expressed  in  terms  of  the  local  surface 
shear.  The  histoiy  of  the  boundary-layer  growth  u~  to  the  point  in  question 
does  not  appear  explicitly,  but  only  through  the  value  of  the  local  shear. 

To  establish  the  position  of  this  local  shear  on  the  surface  of  a  body,  it 
is  necessary  to  first  establish  the  local  boundary- layer  thickness.  This  is 
accomplished  by  linking  the  law  of  the  wall  described  in  the  previous  sec¬ 
tions  with  a  flow  model  applicable  to  the  outer  portion  of  the  boundary  layer. 
This  is  accomplished  in  the  following  sections  of  this  report. 
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2.3  Solution  of  the  Two-Dimens^nnal  Bound  iry-Layer  Equations  (Law  of  the 
Wake) 

As  outlined  in  the  Introduction,  the  turbulent  boundary  layer  is  char¬ 
acterized  by  two  regions;  the  inner  region  that  is  essentially  one  dimensional 
and  the  outer  region  that  is  two  dimensional.  A  general  solution  of  the 
outer,  two-dimensional  region  would  require  an  understanding  of  turbulence 
mechanisms  beyond  the  current  state  of  knowledge  and  an  ingeneous  handling 
of  the  inner  bounc  **y  conditions  to  provide  a  continuous  match  with  the  inner 
region.  Fortunately,  an  empirical  discovery  by  Clauser,  Reference  6,  that 
the  eddy  viscosity  is  essentially  constant  across  the  outer  portion  of  a 
boundary  layer,  but  varies  from  station  to  station,  opens  new  avenues  for  an 
analytical  representation  of  the  wake  region.  In  fact,  as  will  be  shown 
presently,  the  old  standby  in  laminar  boundary -layer  theory  of  local  similar¬ 
ity  can  be  employed  in  this  portion  of  the  turbulent  boundary  layer.  The 

j 

novel  requirements  here  with  respect  to  laminar  boundary-layer  solutions  are 
that  the  inner  boundary  conditions  of  the  two-dimensional  region  must  be 
established  through  a  matching  process  with  the  inner  one-dimensional  law  of 
the  wall  and  that  the  eddy  viscosity  varies  with  tne  distance  along  the  sur- 
face. 

It  will  be  noticed  in  the  following  development  that  the  solution  of  the 
two-dimensional  region  can  be  achieved  through  the  use  of  similarity  argu¬ 
ments.  Because  some  of  the  requirements  for  strict  similarity  are  not  uni¬ 
versally  realized,  it  is  necessary  to  think  in  terms  of  "local"  similarity 
much  as  the  investigators  did  with  laminar  boundary  layers  that  truly  are  not 
similar.  This  mathematical  inexactness  obviates  the  need  for  extreme  thorough¬ 
ness  in  the  matching  of  the  inner  and  outer  flow,  a  piocess  that  would  involve 
iterative  procedures  of  significant  complexity.  Matching  here  is  accomplished 
by  extending  the  two-dimensional  solution  to  the  surface  and  employing  a  fic- 

| 

titious  slip  velocity  there.  This  slip  velocity  becomes  a  variable  that  can 
be  adjusted  in  the  matching  procedure  needed  to  couple  the  flow  solution  for 

I  the  inner  and  outer  solution.  The  matching  procedure  will  be  described  in 
detail  in  Section  2.4.  In  this  section  a  description  is  given  of  the  two- 
dimensional  boundary-layer  solution. 

The  basic  incompressible  boundary-layer  equations  for  which  a  similarity 
solution  is  sought  are  the  continuity  equation.  Equation  (2-6),  reduced  to 

| 

1 
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at  constant  density  and  the  momentum  equation,  Equation  (2-8),  with  the 
inviscid  flow  pressure  gradient  term  retained 


du  ,  du 
u  $ +  v  “ 


7  ai 
m  v  2 
dy 


df» 

+  u 


i  dx 


(2-9:) 


Here 


e  is  the  total  kinematic  viscosity  acting  in  the  fullv  developed 
m 


turbulent  (ox  wake)  region  of  the  boundary  (i.e.,  +  v  has  been  replaced 

by  e  ).  The  form  of  the  third  term  in  Equation  (2-93)  indicates  an  antici- 
m 

pation  of  employing  the  Clauser  criterion  that 


u  6* 

l 


const  ■  0.018 


(2-94) 


namely,  the  dependence  of  on  x  through  6*  and  but  am  independ¬ 

ence  of  the  normal  coordinate  y. 

The  first  step  in  the  formulation  of  a  similarity  solution  involves  a 
selection  of  certain  similarity  forms  for  the  dependent  and  independent 
variables.  For  the  present,  these  are  chosen  as 


u  -  u 


Ui  -  us 


(2-95) 


v  -  v 


v  -  v 
1  o 


(2-96) 


y  .  £ 
y  6 


(2-97) 


where  u  is  the  apparent  wall-slip  velocity  and  is  the  normal  compon- 

ent  of  velocity  existing  at  points  beyond  the  boundary  layer,  that  is,  for 
y  greater  than  6.  In  general,  v  is  a  function  of  y  as  well  as  x,  but 
can  be  related  to  the  other  parameters  of  the  solution  by  use  of  the  continu¬ 
ity  equation,  that  is, 


v  -  v 


jr 

/■ 


(2-98) 
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which  for  y  2  6  yields 


d&*  .  du 

ax  s  lm  dx 


du, 

y  d>r 


(2-99) 


Here  &*sim  is  employed  to  represent  the  displacement  thickness  of  the  entire 
similarity  flow  region  from  y  ■  0  to  y  -  5  and  with  slip  at  the  surface. 

It  should  be  noted  that  by  integrating  Equation  (2-98)  over  the  range  of  y 
down  to  the  wall,  a  region  of  the  boundary  is  being  considered  where  the  outer 
flow  solution  does  not  correspond  to  the  real  flow.  To  account  for  this,  a 
distinction  is  made  between  the  displacement  thickness  evaluated  from  the 

similarity  solution  &*sim  and  the  displacement  thickness  applicable  to  the 
real  flow,  6*. 

Consider  first  the  transformation  of  the  independent  variables  x  and 
y  to  x  and  y. 


(2-100) 


The  transformation  equations  become 


iL  «  _SL  _  £  flfi. 
dx  dx  5  dx  dy 


(2-97) 
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iL  .  1  jL 

dy  5  dy 


(2-102) 


When  the  continuity  equation  is  transformed  according  to  Equation**  (2-101) 
and  (2-102)  and  the  dependent  variables  are  replaced  by  the  definitions 
(2-95)  and  (2-96),  there  results 


(„  -  u  )  fan  _  £  dfi.  A  -  A.  ,  * 

1  s  \d5?  5  dx  dy,/  u  dx  ^ui  ”  us  +  dx 

+  U.  -V,girc  +  (6*  . 

1  dx  su 


i  a.sil 

dx  J  6  ^  dx 


Similar  transformations  convert  the  momentum  equation  to 


(2-103) 
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(2-104) 


If  similarity  exists,  that  is,  G  is  a  function  of  y  only,  the  axial  varia¬ 
tion  of  Q  at  constant  y  must  var.  ih.  Incorporating  this  similarity 
assumption  into  the  shove  relations  and  further  rearranging  their,  results  in 
a  continuity  equation  of  the  form 

-  [k  -  us'  o  if ]?  § +  [i  x  -  = +  +  (-  ^-)v 
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and  a  momentum  equation  of  the  form 
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In  order  to  achieve  similarity,  the  axially  dependent  coefficients  oi  the 
variables  in  the  above  equations  must  be  in  constant  ratio  to  avoid  the 
inconsistency  of  an  x  dependence.  If  the  Clauser  parameter  is  modified  to 
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(a  relative  minor  change)  and  is  presumed  constant,  it  is  required  that  the 
following  parameters  be  axially  invariant  (or  at  least  very  slowly  varying) 

u  v  _  d.u  5*  . 

s  o  d5  5  _ l.  aim 

ux  *  ux  *  dx  *  ut  d5  *  6 


The  first  two  parameters  require  a  linear  dependence  of  the  fictitious  slip 
velocity  and  the  normal  surface  mass  transfer  velocity  on  the  free-stream 
velocity  at  the  edge  of  the  boundary  layer.  For  a  flat  plate,  the  latter 
parameter  would  indicate  the  requirement  of  a  constant  blowing  rate  along 
the  plate.  The  third  parameter  indicates  a  linear  growth  of  the  boundary 
layer.  The  fourth  parameter,  when  combined  with  the  third,  results  in  a 
requirement  on  the  allowed  free-stream  velocity  distribution,  namely - 

ux  *  (ax  +  b)m 


where  a,  b,  and  m  are  constants.  The  fifth  parameter  indicated  ia  not 
really  irdependent  as  can  be  discerned  from  the  definition  of  the  displace¬ 
ment  thickness  expressed  in  terms  of  the  similarity  variables 

-(l  ~Z*)f  <1  -  0)<3y  (2-107) 

o 

The  value  of  the  definite  integral  is  itself  dependent  on  v  /u  and  u  /u.*, 
thus,  the  requirement  of  constant  is  automatically  satisfied  by 

constant  and  u  /uA  and  no  new  requirement  is  imposed. 

The  effect  of  expressing  the  Clauser  parameter  in  terms  of  6* 

slm 

rather  than  5*  is  approximately  5  percent  (i.e.,  e ■  0.018  -* 
€ii/ui6*siin  ■  0.019)  since  the  ratio  of  the  two  displacement  thicknesses  was, 
for  the  entire  range  of  data  studied,  approximately  0.95.  In  the  numerical 
work  reported  here,  this  modification  was  not  employed  and  an  iterative  tech¬ 
nique  was  introduced  to  correct  the  solutions.  It  soon  became  obvious  that 
this  iteration  was  not  warranted  and  that  the  constancy  of  the  modified 
Clauser  parameters  was  as  readily  justifiable  as  was  the  constancy  of  the 
original  parameter. 

*-  ! 
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Because  of  the  requirement  that  u /u,  be  constant,  no  loss  in  general- 

s 

ity  results  from  replacing  Q  with  u/ux  or  u.  Also,  the  requirement  that 

6*  .  /5  be  constant  permits  redefining  y  as  y/&*  These  redefinitions 

sim  r  sim 

result  in  considerable  algebraic  simplification  in  the  transformed  continuity 
equation 


5  *  .  du , 
sim  i 

u.  dx 


(u-v)  - 


d6*  .  /d6*  .  5*  .  du,' 

aim  -  dll  .  ( .  _ fiJJR  .  _ _ L 

dx  y  dy  y  dx  u,  dx 


_ 

dy 


6*  .  y  du. 


u. 


dx 


4 

dy 
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and  in  the  transformed  momentum  equation 


6*  .  du. 
_ SJkBL 


u. 


dx 


v 

U,  dy 


dx 


5*  .  du, 

m  .  _ &2JT  _ 2 

ui  dx 


,•»  d&* 

~  du  ~  sim  ~  du 
dy  1  dx  dy 


6*  . 


u. 


du 

dx 


-f  u2  -  y?  -  — 
V  *  dy  /  ux 


_BL 


d2u 

6*  2 

sim  dy“ 


These  simultaneous  differential  equations,  one  of  the  first  and  the  other  of 
the  second  order,  can  be  combined  into  a  single  third-order  equation.  This 
is  done  by  solving  the  continuity  equation  for  v  and  introducing  it  into 
the  momentum  equation.  Before  doing  this,  however,  it  is  well  to  introduce 
cercain  simplified  nomenclature  by  mear;r  of  the  following  definitions: 


d6< 


a. 


aim 


dx 


u 


u. 


V 

_a 

u. 


pu, 


5*  . 

sim 


u 


a. 


i 

€ 

-JL. 


2L 

dx 


u  5*  . 
l  sim 


2-41 


An  expression  for  the  \f  can  be  obtained  by  integrating  the  continuity  equa* 
tion  once  with  respect  to  u  and  rearranging,  yielding 


(a1  +  a s)  J  y  du  -  a^Q 

________ 


(2-110) 


when  this  relation  is  introduced  into  the  momentum  equation  the  following 
third-order  differential  equation  is  obtained 


d2u 

a  + 

8  —2 

dy 


(ai  +  as)  J  S  dy  -  a4 


df  +  a5U  -  u2) 


(2-111) 


The  integral  in  the  coefficient  of  the  second  term  in  Equation  (2-111)  can 
be  eliminated  by  introducing  a  stream-function  related  term,  namely 


a'  -  ^  -  S 
q  dy  “  u 

In  terms  of  this  variable,  Equation  (2-111)  becomes 

I 

| 

ael‘"  +  [<ai  +  V»  -  »,]  9”  *  V1  -  9'2'  - 

and  the  following  boundary  conditions  apply: 

i 

y 


C-112) 


(2-113) 


/s 


dy  •  0 


(2-114) 


g '  -  u  *  a, 
*  S3 


and 


**  *7 

v  — *  oo  f 


g'  -  u  -  1 


With  the  exception  of  the  term  containing  a4  which  relates  to  surface  mass 
addition,  this  equation  is  identical  to  the  one  presented  by  Ciauser  in 
|  Reference  6. 

Before  the  techniques  employed  in  the  solution  of  the  problem  described 
|  in  Equations  (2-113)  and  (2-114)  are  discussed,  it  is  convenient  to  relate 
the  dependent  variable  g  with  the  actual  shear  at  the  wall  and  with  the 


7See  also  Equation  (2-121) 
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displacement  thickness*  The  momentum  equation  fcr  the  two-dimensional  solu¬ 
tion  with  slip,  Equation  (2-93),  reduces  to  the  following  equation  as  y -*■  0 


Dusl  £u  1  dr  dui 

s  ox  o  oy  oy  i  dx 


(2-115) 


For  u  •  a,u ,  Equation  s2-115)  becomes 
8  3  1 


^  aui  aui  bu 

3y  "  Pa3  ui  d 7  "  dT  +  ?vo 


(2-116) 


which  when  integrated  with  t aspect  to  y  yields 


aui  /  \ 

T  *  Ts  *  PU1  dT  (s2  ~  0  y  +  PVQ(U  "  U8) 


(2-117) 


In  the  actual  boundary  layer,  no  slip  occurs  and  the  corresponding  expression 
for  the  local  shear  is 


t  -  t  ■  pv  u 
o  o 


(2-118) 


For  the  same  shear  condition  away  from  the  surface  in  the  real  boundary  layer 
as  in  the  one  with  slip,  the  following  relationship  is  required  at  y  -  0 


t  +  pv  u  ■  x 
o  r  o  s  s 


(2-119) 


In  terms  of  the  g  function,  Equation  (2-119)  can  be  rewritten  as 


2  pu 


o_.  _  la  + _ (s&\ 

/  ui  ui  ui6*sim 


g"  (0) 


^  +  a3a4 


(2-120) 


Also 


y>o 

■  f  <l  -  3)dy 


*Note  that  the  additional  local  similarity  req.ij.rement  that  Cf/2  be  constant 
is  introduced  here. 


li  iih  writ 
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-/  '"s 


-  u)dy  ■  y  -  g 


g  -  y  -  1 


(2-121) 


at  large  y.  Equation  (2-121)  represents  an  additional  boundary  condition 
for  large  y  in  Equation  (2-114) . 

To  avoid  the  iterations  inherent  in  the  numerical  integration  of  the 
boundary-value  problem  discussed  by  Equations  (2-113)  and  (2-114),  the  basic 
problem  is  transformed  to  an  initial  value  problem  through  the  transformations 


f  -  Clg  ^ 

i  m  cay  j 


(2-122) 


where  cx  and  c2  are  constants.  Restricting  Equation  (2-113)  to  no  main 
stream  acceleration,  that  is,  a  flat  plate  with  as  -  0,  the  transformations 
(2-122)  result  in 


Af*  -  BffH  -  f •  1 


(2-123) 


with  the  boundary  conditions 


(2-124) 


where  the  parameters  A,  B,  C,  D,  E,  and  the  function  F  are  defined  by 


c.a,. 

2  fi 


(2-125) 


a.  +  as 

_ 34 


/  f 

c,  !  ~ 
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It  can  bo  seen  that  if  an  arbitrary  set  A,  B,  C,  and  D  is  specified  then 
a  complete  initial  value  problem  hr.s  been  formulated.  The  solution  of  this 
problem  v/ill  yield  the  parameter  E  and,  from  the  def iritions,  the  function 
F.  It  is  seen  that  knowledge  of  these  latter  two  quantities  permits  the 


evaluation  of  c. 


and 


and  subsequently  the  evaluation  of  the  remaining 


'i  -  *2 

parameters.  From  this  process  there  results  a  self  consistent  set  of  values 
for  the  parameters  a  ,  a4/a8,  a  /a  ,  and  a2/a0  without  recourse  to  iterat¬ 
ing  the  solution  of  the  differential  equation.  Thus,  in  terms  of  physical 
quantities  a  consistent  set  is  found  between  C^/2,  v Q/ui,  u  /u^  and 

d6Wdx- 

Equation  (2-123)  was  solved  using  a  standard  Runge-Kutta  routine.  To 
test  this  parametric  procedure,  it  was  applied  to  the  Blasius  problem.  For 
this  case  of  a  laminar  boundary  layer  aa  is  no  longer  a  constant*,  however,9 


and 


ai  d  Re&* 

T  -  const  "  Re&*  dR e  “  c2A 

e  x 


ai  5f  c2‘ 

—  »  const  -  Re  -  *  ■ 

cl  /  O  C . 

O  JL 


(2-126) 
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Also,  with  zero  mass  addition  and  surface  slip,  there  results 
Thus,  the  problem  is  defined  as  an  initial  value  problem  when 
are  arbitrarily  selected.  Two  independent  sets  of 
ized  and  both  yielded 


3 

A,  b, 


a 


4 

and 


0. 


A,  B,  and 


It 


0.33205  Re 


-1/2 


C  were  util- 


(2-128) 


the  classical  laminar  boundary-layer  solution. 

In  the  numerical  work  associated  with  the  similarity  portion  of  the 
boundary-layer,  values  of  A,  B,  C,  and  D  were  arbitrarily  assigned.  I', 
was  found  that  72  independent  solutions  were  adequate  to  provide  solutions 
covering  the  following  ranges: 


}If  after  Equation  (2-106),  €m  is  considered  axially  independent,  the  result¬ 


ant  set  of  ax’^lly  independent  parameters  is  a  and  u.o1 


a2»  a4> 


and 


sim 


times 


li » 


e  * 
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0.3  <  us/ux  \  1.0 
-0.003  v  vQ/ux  v  0.005 

and  a  wide  variety  of  values  of  Cf/2  and  d5*/dx.  How  these  results  are 
related  to  the  "law-of-the-wall"  regio’  will  be  described  in  detail  in  the 
next  section. 

j  2.4  Matching  Procedure  for  Solutions  of  the  Laws  of  the  Wall  and  the  Wake 

In  the  previous  sections,  solutions  of  the  laws  of  the  wall  and  wake 
were  found  having  the  following  functional  relationships 

Law  of  the  wall  u+  ■  F(y+,  vq+)  (2-129) 

|  Law  of  the  wake  u  ■  G(y,  Vg/u^  Cf/2,  Ug/u.^  do*/dx)  (2-130) 

The  four  parametric  terms  on  the  right  of  Equation  (21-30)  (other  than  y) 
are  not  randomly  variable  but  were  found  to  be  grouped  in  self  consistent 
sets  in  which  two  of  the  parameters  may  be  considered  independent. 

Several  matching  procedures  are  possible  to  3  ink  these  solutions,  however, 
the  one  described  here  is  thought  to  be  the  most  straightforward  in  concept. 

The  matching  procedure  is  begun  by  selecting  arbitrarily  one  of  the 
sets  of  parameters 

I  In  2*.  flSl 

f  u  *  2  *  u  *  dx 

%  i  i 


found  to  be  self  cox  it.  The  profile  results  of  the  two-dimensional 

similarity  solution  i^  .»en  plotted  as  u  versus  y  “  y/5*  Knowledge 

of  C^/2  and  permits  the  direct  evaluation  of  v^*  from  Equations 

(2-33)  and  (2-36)  as 


uiVv2 


(2-131) 


Identification  of  \ 
ship  of  u+  versus 


permits  establishing  the  particular  unique  relation¬ 


ship  of  u  versus  y  corresponding  to  Equation  (2-130)  for  the  inner 
region.  At  this  point  a  particular  'Snatch  and,  thus,  u*atch  is  tentatively 
chosen  as  the  match  point.  Since  the  value  of  e ^/v  increases  continuously 


tM««  . . . . . . . WWW 
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match 


aiso 


from  unity  at  y  *  0  all  along  the  inner  solution,  chosing  y 
yields  a  value  of  (~n/v)mat.ch*  s*nce  the  outer  region  was  assumed  to  have 
a  constant  value  of  T  ,  it  is  the  value  of  ^”r/v^match  fc^at  18  applicable 
over  the  entire  outer  region.  From  Equation  (2-94)  it  can  be  shown  that 


Re 


b*. 

sim 


u  b* . 

1  sun 

v 


(V'O. 


match 


(2-132) 


Utilizing  this  information,  it  is  easy  to  relate  the  coordinates  of  the  inner 
and  outer  flows,  namely. 


u  -  u+yCf/2 


and 


V 


Re&*  VV7 
sxm 


(2-132) 


(2-134) 


The  value  of  ymatch  is  moved  about  iteratively  until  u^^  and  ymatch 
lie  on  the  u(y)  curve  of  the  outer  flow  solution.10 

The  resulting  curve  of  u(y)  is  then  utilized  to  evaluate  local 
boundary -layer  thickness  characteristics  needed  in  subsequent  correlations. 
The  displacement  thickness  of  the  confined  boundary  layer  is  given  by 


Xmatch 


(2-135) 


and  the  momentum  thickness  is  given  by 


^match 


u(l  - 


sim 


-  u)  dv 


(2-136) 


10Although  the  matching  method  is  described  as  a  graphical  process,  in  actual 
practice  the  procedure  was  programed  to  be  performed  by  the  computer.  A 
mismatch  in  the  local  shear  also  results;  however,  this  has  been  estimated 
as  less  than  0.5  percent  of  the  local  shear. 
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The  integrals  utilize  the  velocity  distribution  appropriate  to  the  inner  and 
outer  layers  of  the  boundary  layer.  Equations  (2-135)  and  (2-136)  combined 
with  Equation  (2-132)  yield  the  displacement  and  momentum  thickness  boundary 
layer  appropriate  to  the  blowing  parameter,  Vg/u^  and  skin-friction  coef¬ 
ficient  C^/2,  in  the  Equation  (2-131). 


2.5  Evaluation  of  Reynolds  Numbers  Based  on  the  Length  of  Run 


§ 

| 

I 


I 


The  previous  section  described  the  technique  for  establishing  the 
relationships 


(2-137) 


(2-138) 


for  the  entire  boundary  layer.  To  relate  the  skir.-friction  coefficient  with 
the  Reynolds  number  based  on  the  length  of  run,  use  is  made  of  the  integral 
form  of  the  momentum  Equation  (2-51).  For  a  flat  plate,  this  equation  becomes 


v 

-Q. 


(2-139) 


For  constant  vQ/ui  and  turbulent  flow  from  the  leading  edge,  this  equation 
can  be  integrated  conveniently  to  yield 


(2-140) 


The  integral  in  Equation  (2-140)  was  evaluated  numerically  from  the  results 
indicated  by  Equation  (2-137).  In  the  vicinity  of  the  leading  edge  it  was 
necessary  to  extrapolate  to  Rq  smaller  than  given  by  the  numerical  examples 
worked  out.  This  was  accomplished  by  first  deriving  an  exponential  equation 
for  Cf/2  in  terms  of  Re^  from  the  small  Rey  numerical  results  and, 
thus,  using  this  ecuation  for  extrapolating  the  very  small  values  of  R  . 

If  one  accepts  the  argument  that  Equations  (2-137)  and  (2-138)  are  the 
basis  of  a  local  similarity  even  under  conditions  where  streamwise  pressure 
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gradients  occur,  the  entire  momentum  equation,  Equation  (2-51),  can  be  solved 
numerically  to  also  yield  the  growth  of  the  boundary  layer  on  two-dimensional 
bodies  other  than  flat  plates.  It  would  be  preferable,  however,  to  expand 
the  set  of  independent  parameters  to  include  the  local  pressure  gradient 
parameter,  a_.,  and  then  enter  the  integral  momentum  equation. 

2.6  Solution  of  the  Energy  and  Diffurion  Equations 

The  solutions  of  the  energy  and  diffusion  equations  follow  nearly  the 
same  course  as  did  the  solution  of  the  momentum  equations,  except  that  the 
matching  of  the  innner  and  outer  solutions  is  accomplished  at  a  known  posi¬ 
tion  in  the  boundary  layer  where  momentum  matching  was  already  established. 
The  restriction  to  the  constant  property  case  in  the  current  analysis  permits 
separating  the  solutions  of  the  individual  equations  and  utilizing  the  exist¬ 
ing  momentum  solution  to  solve  the  energy  and  diffusion  equations.  In  the 
more  general  case,  these  scalar  ccnservational  equations  would  have  to  be 
coupled  with  the  momentum  equation  and  simultaneous  solution  of  the  three 
equations  would  be  required. 

The  basic  energy  conservation  equation  for  a  constant  property  boundary 
layer  is 


u  aX  .  v  ax  .  _m_  afx 

Z*  ^  Pr  ’  dy2 


(2-141) 


where  the  Prandtl  number,  Pr',  is  defined  here  as  the  overall  ratio  of  momen¬ 
tum  to  energy  transfer  coefficients,  that  is. 


Pr'  - 


e_  +  v  c 

_ D _  »  m 

e  +  JL.  €  +  Js_ 

pC  pC 

P  P 


(2-142) 


again  replacing  the  sum  of  €m  +  v 


The  value  of  Pr'  in  the  two- 


dimensional  region  is  obtained  from  the  one-dimensional  energy  solution  at 

the  y  +  determined  in  the  matching  of  the  one-  and  two-dimensional  momentum 
m 

solutions.  Generally,  at  the  match  point  the  contribution  of  the  molecular 
mechanisms  is  negligible  compared  vo  the  turbulent  processes  and 


~  Id 


(2-143) 
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Continuity  of  €ffl,  Pr ' ,  temperature,  and  temperature  derivatives  across  the 
match  point  assures  a  continuity  in  the  heat  flux  at  that  point.  For  the 
diffusion  equations,  all  relations  would  be  identical  with  replacing  T 

and  Sc'  replacing  Pr'.  Here, 


e  +  v 

Sc'  -  -m - 

+  n 
D 


(2-144) 


When  the  outer  flow  is  considered,  it  is  convenient  to  define  the  depend¬ 
ent  variable  os 


_  T  -  T 

T  *  T  T  "  1  -  T 
O  1 


(2-145) 


When  Equation  (2-145)  is  inserted  into  Equation  (2-141)  and  the  transformations. 
Equations  (2-101)  and  (2-102),  are  utilized,  there  results 


where 


“a 
Pr ' 


iOL 

a?2 


[(*!  +  -  “,] 


-  a  uT  *  0 
dy  7 


(2-146) 


a 
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(2-147) 


Similar  expressions  for  the  diffusion  equation  are 


where 


,~2 

dy 


[(ax  +  a5)g 


d? 


a  uX 
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(2-148) 


dx 


(2-149) 


The  boundary  conditions  on  these  equations  are  the  requirement  of  a  match 
with  the  inner  region  at  ymatch  and  that  T  -*  0  and  X  -*  0  as  y  — 

In  both  of  the  above  equations,  it  was  assumed  the  free-strearo  properties, 
T!  and  ar®  invariant  with  distance  along  the  surface  of  the  body. 

For  consistency  with  the  previous  momentum  equation  solutions,  a  must  be 
set  equal  to  zero. 
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For  the  case  of  a  constant  surface  temperature,  a^  ■  0,  and  Equation 
(2-146)  reduces  to 


+  <^9  *  a,>  if  *  0 


(2-150) 


When  Equation  (2-150)  is  integrated  with  respect  to  y,  there  results 
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atch 


which,  when  integrated  again,  yields 


T  -  T 
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dy  (2-152) 


^match 


At  the  upper  boundary  condition,  there  results 


-  f  (ai9-g4)dy 


match 


^match 
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dy  -  I 


(2-153) 


'match  ^match 


Since  g(y)  is  known  from  the  solution  of  the  momentum  equation,  the  integral 
in  Equation  (2-153)  can  be  evaluated  directly  fiom  the  chosen  self-consistent 
set  of  d&*/dx,  Cj/2,  u^/u^  v Q/u1,  and  the  known  match  point. 

From  the  solution  of  the  inner  region  T+(y+,  vq+) ,  and  the  definition, 
Equation  (2-5a),  it  is  seen  that  that 


T  )  - 


Sl_  t+ 


Vv^ 


(2-154) 


and 
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m  -  r  R  * 

dy  Si  o* 


. 

0*  .  ,  + 

sim  dy 


(2-155) 


From  a  combination  of  Equations  (2—153)  through  (2-155),  it  is  found  tnai 


(2-156) 


match 


match 


All  the  terms  of  the  right  member  of  Equation  (2-156)  are  readily  found  from 
the  prescribed  values  of  Cf/2  and  v  and  ymatc^  established  in  the 

momentum  solution.  From  Equations  (2-154)  through  (2-156)  evaluation  of 
Thatch  and  (dT/dy)  atc^  can  be  ma<*e  anc^  ***th  these  quantities,  the  tem¬ 
perature  distribution  through  the  boundary  layer  is  given  by  Equation  (2-152) 
and  the  numerical  inner  flow  solution  T  (y  ) .  The  energy  thickness  (energy 
decrement  thickness)  is  found  by  numerically  integrating  the  following 


match 


uT  dy  ♦ 


uT  dy 


(2-157) 


0 

*  match 


over  the  two  regions  of  the  boundary  layer.  The  length  Reynolds  number 
corresponding  to  Cjj  and  £  found  in  this  manner  is  the  same  corresponding 
to  the  Cf/2  used  as  a  parameter  in  the  above  procedure. 

For  the  case  of  constant  v along  a  plate,  the  most  common  injec¬ 
tion  distribution  utilized  in  the  existing  experiments,  a  chtnge  in  xq,  the 
molar  concentration  at  the  surface,  occurs  with  distance  along  the  surface. 
This  change  in  xq  results  in  the  requirement  for  a  non-zero  value  of  afl 
(Eq.  (2-149))  in  Equation  (2-148).  Under  this  condition,  Equation  (2-148) 
cannot  be  solved  in  closed  form  and  numerical  integration  is  required.  To 
avoid  iterations  in  the  integrations,  advantage  is  taken  of  the  homogeneous 
nature  of  Equation  (2-148)  to  express 


where  at  large  y 
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Equation  (2-148)  thus  becomes 
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which  can  be  integrated  for  specified  values  of  Cf/2  and  and  assumed 

values  of  a6.  The  boundary  conditions  that  muse  be  satisfied  at  ymatch 
(continuity  of  concentration  and  i.tr  derivative)  require  that 
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The  value  o'  CM  is  found  cor-’eniently  from 
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Both  a  and  C^,  expressed  above,  are  dependent  on  the  choice  of  a0  util¬ 
ized  in  Equation  (2-159).  The  unique  value  of  a0  appropriate  to  the  solu¬ 
tion  is  found  iteratively  from  the  requirement  that  over  a  plate  with  uniform 
blowing 
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is  evaluated  numerically.  The  parameter  a@  is  iterated  until  the  Rep 
found  from  Equations  (2-162)  and  (2-163)  agree. 

A  similar  situation  to  the  above  occurs  for  the  energy  boundary  layer 
on  a  porous  plate  dissipating  a  constant  heat  flux  with  a  uniform  coolant 
gas  injection.  Here  a^  is  no  longer  zero  and  has  to  be  iterated  so  that 
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(2-164) 


where  q  is  the  net  energy  supplied  to  the  wall,  for  example,  electrical 
heating  less  radiation  losses  plus  bulk  convection  ‘pv0cp(T2  “  Tx))* 
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3.  RESULTS  AND  D7SCUSSI0F 

In  the  preceding  section  a  set  of  analytical  procedures  was  developed 
for  the  prediction  of  the  profiles  of  velocity,  temperature,  and  concentra¬ 
tion  in  a  turbulent  boundary  layer.  The  results  predicted  by  the  " law-of- 
the-wall"  relations  were  compared  with  data  in  order  t:o  establish  the  valid¬ 
ity  of  that  portion  of  the  model.  For  the  outer  (or  wake)  region  of  the 
boundary  layer  a  locally  similar  solution  was  devised.  When  this  latter 
solution  is  coupled  with  the  inner  "wall  law"  solution  the  complete  profile 
predictions  result.  These  profiles  can  subsequently  serve  as  the  basis  for 
integral  solutions  of  the  boundary  layer. 

In  this  section  the  pr  'dictions  for  the  complete  profiles  will  be  com¬ 
pared  with  data  and  the  subsequent  results  of  the  integr?’  solutions  will  be 
presented,  compared  with  data,  and  discussed. 

However,  one  new  postulate  was  introduced  in  the  .nalysis  of  the  outer 
region  of  the  profile  that  should  be  evaluated  prior  to  the  presentation  -»f 
the  results  obtained.  This  postulate  is  related  to  the  Clauser  parameter  and 
its  constancy  under  conditions  of  surface  mass  transfer. 

The  experimental  variation  of  eddy  viscosity  across  a  boundary  layer 
where  significant  surface  blowing  occurs  is  shown  in  Figure  25  along  with 
hot-wire  data  taken  over  a  smooth,  solid  plate  (the  filled-in  symbols).  The 
eddy  viscosity  corresponding  to  data  with  blowing  was  evaluated  by  employing 
Equation  (2-48)  and  velocity  profiles  from  Reference  17.  Although  this  latter 
technique  was  not  expected  to  be  very  accurate,  rather  good  agreement  was 
obtained  with  the  more  accurate  hot-wire  data.  The  data  for  the  blowing  case 
exhibit  a  scatter  that  lies  almost  equally  above  and  below  the  more  accurate 
data.  A  value  of  the  Clauser  parameter,  eR/ux^*»  of  0.018  agrees  quite  well 
with  both  sets  of  data  for  0.2  <  y/5  <  0.6.  Above  y/5  *  0.6  the  parameter 
drops  oft,  primarily  due  to  intermittency ,  that  is,  the  fraction  of  the  time 
the  boundary  layer  is  turbulent  at  a  particular  value  of  y/5.  That  this  is 
the  case  is  seen  by  the  "corrected"  data  shown  on  Figure  25  which  were  obtained 
by  dividing  the  eddy  viscosity  values  by  the  values  of  the  intermittency 
factors  given  in  Reference  51.  The  most  .significant  conclusion  that  can  be 
drawn  from  this  figure  i<*  chat  e^/u^*  can  be  considered  to  be  independent 
of  surface  blowing  and  this  is  the  basis  of  Equation  (2-94). 

A  comparison  is  made  in  Figure  26  between  experimental  velocity  profiler 
measured  over  smooth,  non-porous  plates,  and  the  predicted  results  of  the 
current  analysis.  The  theoretical  praaictions  nhown  are  evaluated  at  che 
value  of  Re^  that  corresponded  to  the  data.  This  basis  of  comparison  on 
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the  local  boundary-layer  thickness  avoided  those  uncertainties  that  occur 
when  the  length-of-run  Reynolds  number  is  employed,  namely,  accounting  for 
the  effects  of  the  position  of  natural  transition  downstream  of  the  plate 
leading  edge  or  of  boundary-layer  trips.  In  Figure  26(a)  the  velocity  pro¬ 
file  comparison  is  made  with  the  data  of  Reference  51.  Except  for  the  data 
points  quite  close  to  the  surface,  the  prediction  for  velocity  at  a  given 
distance  from  the  surface  is  within  a  few  percent  of  the  data.  In  the  inner 
region,  the  difference  between  the  prediction  and  the  data  becomes  as  large 
as  10  percent  where  the  absolute  values  of  the  velocity  ar'-  much  smaller.  No 
systematic  difference  between  the  analytical  results  and  the  data  can  really 
be  expected  because  the  Rep  were  matched*,  however,  the  prediction  over¬ 
predicts  and  under-predicts  in  four  alternate  regions  indicating  a  generally 
good  agreement  with  the  data  over  the  whole  of  the  boundary  layer  rather  than 
favoring  one  boundtry-layer  region  over  the  other.  The  comparison  shown  in 
Figures  26(b)  and  (c)  with  the  data  of  References  52  and  11  indicates  essen¬ 
tially  an  identical  behavior. 

In  Figure  27,  similar  velocity  profile  comparisons  are  made  with  data 
obtained  over  porous  plates  where  surface  blowing  or  suction  occurred.  Again, 
the  theoretical  and  experimental  values  of  Re&  were  matched.  In  Figure  27(a) 
comparison  of  the  analytical  results  is  made  with  the  data  of  References  16 
and  17,  which  data  cover  a  large  range  of  uniform  blowing  and  suction  condi¬ 
tions.  Again,  the  general  agreement  between  the  analysis  and  the  blowing 
data  is  quite  good,  the  differences  in  velocity  being  of  the  order  of  a  few 
percent,  except  for  the  inner  2  percent  of  tne  boundary-layer  thickness  where 
the  prediction  is  about  10  percent  lower  than  the  data.  The  significance  of 
the  latter  observation  is  tempered  when  it  is  remembered  that  the  probe  size 
was  rather  large  compared  to  the  distance  away  from  the  surface  in  this  region 
so  that  data  errors  can  be  expected  to  be  larger  here.  Tht  comparison  with 
the  suction  data,  the  data  points  indicated  by  square  symbols,  is  a] so  very 
good.  Except  for  one  data  point,  differences  of  less  than  4  percent  occur 
over  the  entire  boundary  layer.  The  data  of  Reference  11,  shown  in  Figure 
27(b),  show  the  same  general  agreement  with  the  predictions  as  do  the  data 
of  Reference  17  discussed  previously.  In  the  immediate  vicinity  of  th?  wall, 


*• 

the  differences  between  the  data  and  predictions  become  quite  large  and,  again, 
probe  size  effects  could  be  important  here  as  well  as  the  dynamic  pressure 
contributions  of  turbulence  discussed  in  Section  2. 2. 3. 4. 

It  can  be  concluded  from  Figures  26  and  27  that  the  analysis  does  quite 
a  good  job  in  predicting  the  shapes  of  velocity  profiles  under  conditions  of 

**  'I 

zero  mass  transfer,  blowing,  and  suction. 
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A  comparison  of  the  analytical  prediction  of  temperature  distributions 
with  data  is  presented  in  Figure  28.  These  data  taken  from  Reference  12, 
apply  for  the  case  of  zero  blowing  and  for  ■  0.004  and  the  predic¬ 
tions  are  based  on  matching  the  analytical  with  the  local  measured  values  of 

Re^ .  The  agreement  between  the  predictions  and  the  data  is  within  3  percent 

over  the  entire  boundary,  except  for  the  innermost  points.  From  this  figure 
it  can  be  concluded  that  the  proposed  analytical  model  predict?  these  limited 
number  of  cases  of  temperature  distribution  data  quite  well  under  the  condi¬ 
tions  of  zero  mass  transfer  or  with  blowing. 

The  distribution  of  helium  that  was  injected  into  a  boundary  layer, 
Reference  17,  is  compared  in  Figure  29  with  the  prediction  given  by  the  pres¬ 
ent  analysis.  Again,  the  measured  Re^  is  used  as  the  matching  criterion 

for  the  analysis.  Data  are  presented  for  a  case  with  very  little  blowing, 
v (J/ui  *  0.000113,  and  one  with  rather  large  blowing,  v  /u1  *  0.00545.  For 

the  low  blowing  case,  the  agreement  between  the  analysis  and  the  experiment 

is  within  5  percent,  except  for  the  data  points  right  next  to  the  surface. 
Again,  probe  size  may  have  an  influence  here.  At  the  higher  blowing  rate, 
the  data  are  generally  lower  than  the  predicted  profile,  with  a  maximum  dif¬ 
ference  of  about  10  percent.  Generally,  the  differences  are  much  lower  than 
this  maximum  figure.  It  can  be  concluded,  therefore,  that  the  proposed  tech¬ 
nique  does  a  satisfactory  job  in  predicting  concentration  distribution*. 

The  analytical  results  for  the  model  of  turbuler**  boundary  layer  utilized 
in  this  analysis  are  presented  in  terms  of  dimensionless  transfer  coefficients 
and  shape  factors  as  functions  of  Reynolds  numbers  and  transpiration  rates 
in  Figures  30  through  36. 

Figure  30  shows  the  local  skin-friction  coefficient  on  a  flat  plate  with 
uniform  injection  as  a  function  of  the  Reynolds  number  based  on  the  momentum 
thickness.  The  transpiration  rate,  is  shown  as  a  parameter  of  thc- 

individual  curves.  For  the  curves  representing  suction,  negative  a 

limiting  line  is  shown  where  the  local  suction  rate  is  just  equal  to  the 
local  skin-friction  coefficient,  which  then  results  in  a  constant  boundary- 
layer  thickness  beyond  that  point  (line  of  asymptotic  suctior  solutions). 

The  curves  shown  of  this  figure  reveal  that  outward  transpiration,  positive 
vq/u, ,  reduces  the  local  skin  friction  Rey  relationship,  with  the  reduction 
being  greater  at  the  larger  values  of  Rew •  Suction  increases  the  skin- 
friction  coefficient.  These  results  are  qualitatively  compatible  with  past 
findings . 
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In  Figure  30(b)  the  v  /u  ■  0  line  from  Figure  30(a)  is  reproduced 
along  with  the  recent  data  from  Smith  and  Walker,  Reference  50.  The  maximum 
difference  between  the  data  and  the  prediction  is  5  percent  over  the  Reynolds 
number  range  3xl03  to  5xl04  while  the  bulk  of  these  data  agree  with  the 
prediction  to  a  few  percent. 

Tie  theoretically  predicted  variation  of  the  boundary-layer  form  factor, 

H  «  6V^»  is  shown  in  Figure  31(a)  as  a  function  of  the  momentum-thickness 
Reynolds  numbers  at  a  variety  of  transpiration  rates.  The  curves  on  the 
figure  reveal  that  H  varies  quite  markedly  with  both  Reynolds  number  and 
transpiration  rate.  Consideration  of  this  variation  is  of  considerable 
importance  in  the  evaluation  of  the  turbulent  boundary- layer  growth  on  a  sur¬ 
face  possessing  a  streamwise  pressure  gradient,  see  Equation  (2-51). 

In  Figure  31(b)  the  predicted  variation  of  H  under  zero  transpiration 
is  compared  with  the  data  of  Reference  50.  This  comparison  indicates  that 
the  boundary-layer  model  proposed  herein  satisfactorily  represents  the  form 
factor  under  conditions  of  zero  transpiration  to  within  about  4  percent. 

The  inadequacy  of  power  law  formulations  that  predict  constant  values  of  H 
with  Reynolds  number  is  evident  here. 

When  the  von  KArmAn  integral  equation  is  evaluated  on  a  flat  plate  util¬ 
izing  the  Cf/2(Re0)  relationships  of  Figure  30,  the  curves  of  C^/2 (Re^) 
shown  in  Figure  32  result  for  a  boundary-layer  turbulent  from  the  leading 
edge.  The  values  of  Cf/2 (R e^)  utilized  at  very  low  Re^  were  obtained 
from  logarithmic  extrapolations  of  the  curves  of  Figure  30(a).  It  should  be 
noted  that  the  v  /ux  *  0  line  is  within  a  few  percent  of  tne  KArmAn- 
Schoenherr  relationship  over  the  entire  range  of  Reynolds  numbers  shown. 

Again,  the  increased  effectiveness  of  uniform  t..  <;nrpiratirn  at  the  higher 
Reynolds  numbers  is  evident. 

Similar  curves  to  those  described  above  for  the  heat-transfer  coefficient 
(Stanton  number)  and  the  mass-transfer  coefficient  with  uniform  transpiration 
on  a  flat  plate  are  presented  in  Figures  33  to  36.  The  heat-transfer  coef¬ 
ficients  apply  for  a  laminar  Prandtl  number  of  0.71  and  a  turbulent  Prandtl 
number  of  0.75.  The  mass-transfer  coefficients  apply  for  laminar  and  turbu¬ 
lent  Schmidt  numbers  of  0.2  and  0.75,  respectively.  Again,  the  effect  of 
increared  transpiration  is  greater  at  the  higher  Reynolds  numbers. 

Anothei  way  of  demonstrat ’ ng  the  effect  of  surface  mass  transfer  ori  the 
local  skin-friction  coefficient  is  shown  in  Figure  37.  In  this  coordinate 
system,  C^/C^*  is  plotted  against  v o/u1(C^/2)*  where  the  (*)  refers  to  no- 
blowing  conditions  at  the  same  Reynolds  number.  The  effect  of  !i  'th  Reynolds 
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number  is  largely  eliminated  on  such  a  plot,  as  is  indicated  in  the  current 
results  and  those  of  Reference  22.  The  current  result.*,  are  essentially  the 
same  as  those  of  the  pioneer  paper  by  Dorrance  and  Dore ,  Reference  21,  and 
slightly  lower  than  the  results  of  Rubesin,  Reference  22. 

There  has  been  some  argument  as  to  the  appropriateness  of  the  coordinate 
system  employed  in  Figure  37  because  the  term  v^u^ C^/2*  is  more  consist¬ 
ent  with  the  law  of  the  wall,  Reference  54.  This  parameter  was  used  as  the 
abscissa  of  the  curves  shown  in  Figure  38  and  a  decided  Reynolds  number  effect 
is  introduced.  From  a  correlation  viewpoint,  therefore,  the  coordinates  of 
Figure  37  are  more  appropriate.  Apparently  the  wake  legion  has  a  significant 
effect  in  the  Reynolds  number  and  blowing  parameter  range  shown. 

In  the  evaluation  of  skin  friction  on  bodies  other  than  a  flat  pla' ., 
the  effect  of  surface  mass  transfer  at  constant  values  of  Re^  is  required 
and  this  is  shown  in  Figure  39.  Again,  the  choice  of  the  coordinates  elimin¬ 
ate  the  large  Reynolds  number  dependence.  At  a  given  value  of  the  blowing 
parameter,  the  reduction  in  skin  friction  shown  in  this  figure  is  less  than 
that  shown  in  Figure  37.  This  is  consistent  for  a  boundary  layer  whose  thick¬ 
ness  grows  at  a  rate  slightly  less  than  linear  with  distance  along  a  plate. 

Figure  40  shows  the  effect  of  surface  mass  transfer  on  the  local  mass- 
transfer  coefficient.  The  curves  are  generally  similar  to  the  skin-friction 
results  shown  in  Figures  37  and  39.  Again,  similar  results  for  the  local 
heat-transfer  coefficient  are  given  in  Figure  41. 

The  effect  of  surface  mass  transfer  on  the  coefficients  of  diffusi  n 
and  heat  transfer  can  best  be  shown  by  the  Reynolds  analogy  factors  Cc/2CU 
(dashed  lines)  and  (solid  lines;  indicated  in  Figure  42.  Whereas 

in  Figure  37  there  was  little  to  be  gained  by  using  the  current  analysi.s  in 
place  of  the  early  heuristic  theories,  in  this  figure  it  is  noted  that  the 
heat  transfer  at  the  highest  blowing  rate  is  only  0.4  that  given  by  Refer¬ 
ence  21  (where  Pr  *  1,  Prt  ■  1)  and  0.5  of  that  given  by  Reference  22  (where 
Pr  ■  0.75,  Prfc  *  1).  These  major  differences  result  largely  from  the  more 
precise  boundary-layer  profiles  employed  in  the  current  analysis.  Of  interest 
also  is  that  at  zero  mass  transfer  the  Reynolds  analogy  for  heat  transfer 
agrees  quite  well  with  the  old  Colburn  value  of  Pr2^3.  As  surface  mass 
tran-i-ar  increases  the  Reynolds  analogy  factor  drops  significantly,  indicat¬ 
ing  a  les^r  effect  of  blowing  on  heat  transfer  than  on  skin  friction. 

Because  of  the  very  low  Schmidt  number  (Sc  *  0.21)  for  helium  at  small  con¬ 
centrations,  diffusing  in  air,  the  Reynolds  analogy  factor  for  diffusion  is 
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much  lower  than  that  for  heat  transfer,  although  not  as  low  as  Sc2//j 
zero  blowing.  Again,  surface  mass  transfer  reduces  this  factor  significantly, 
indicating  the  lesser  sensitivity  of  the  mass  transfer  or  diffusion  coeffi¬ 
cient  to  surface  mass  transfer.  For  the  diffusion  coefficients,  the  coordi¬ 
nates  employed  in  Figure  42  do  not  completely  eliminate  the  Reynolds  number 
effect. 


4 .  CONCLUDING  REMARKS 


The  results  of  this  report  demonstrate  that  a  model  of  the  turbulent 
boundary  layer  with  surface  mass  addition  car.  be  established,  that  is,  in 
excellent  agreement  with  profile  data  for  velocity,  temperature,  and  concen¬ 
tration  in  constant  property  boundary  layers  (low  speed,  nmall  temperature 
and  concentration  differences  between  surface  and  free  stream) .  A  uniquely 
interesting  result  of  this  analysis  is  that  the  four  parameters  employed, 
namely:  (1)  the  Prandtl  mixing  length  parameter  K;  (2)  the  Clauser  eddy 

diffusivity  Reynolds  number;  (3)  the  inner  mixing  length  parameter  y  +, 
analogous  to  Rotta's  intercept  parameter;  and  (4)  the  turbulent  Prandtl  or 
Schmidt  number,  can  be  treated  as  constants  independent  of  the  rate  of  sur¬ 
face  mass  transfer  and  Reynolds  number. 

This  agreement  with  the  data,  restricted  though  they  may  be  to  essen¬ 
tially  constant  properties,  is  sufficiently  encouraging  to  suggest  the  next 
step  in  this  work.  Velocity  profile  data  have  been  obtained  in  boundary 
layers  in  supersonic  flow  with  inert  gas  mass  transfer  at  the  surface,  Refer 
enas  5  6  and  55.  It  would  be  a  logical  next  step  in  this  research  to  deter¬ 
mine  the  modifications  needed  in  the  present  theory  to  permit  its  extension 
to  include  the  effects  of  variable  fluid  properties  anu  „o  yield  agreement 
with  the  existing  data.  In  accomplishing  this  extension,  it  will  be  worth¬ 
while  to  examine  the  applicability  of  several  of  the  provocative  ideas  of 
Reference  57  toward  solving  the  current  problem.  Another  phase  of  any 
extensions  of  this  work  must  be  the  planning  of  definitive  experiments  where 
chemical  reactions  occur  as  a  consequence  of  the  surface  mass  transfer.  Pro 
file  data  do  not  exist  at  present  for  such  conditions  and  it  is  only  upon 
such  datt.  that  sound  extensions  of  this  work  can  be  based. 
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(a)  Velocity  distribution  in  a  pipe,  Reference  48,  Reynolds  number  based  on  diameter  -  5.0xl0b 

Figure  4--  Comparison  of  mixing  length  formulations  with  data. 


I 


I 

I 

1 

I 

I 

I 


m  <N  csj 


n 

+ 


(b)  Velocity  distribution  over  a  flat  plate,  Reference  49. 

Figure  4.-  Contir  led. 
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Figure  4.-  Concluded. 
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Figure  7.-  Comparison  of  predicted  direct  viscous  dissipation  and  turbulent  energy 
production  rates  with  the  data  of  Klebanoff  and  Laufer. 


Figure  8.-  Variation  of  sublayer  thickness  with  transpiration  rate  for  two 
hypotheses  in  the  two-layer  boundary-layer  model. 
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figure  9.-  fredictsd  variation  of  di'^wasionless  velocity  with  and  without  transpiration. 
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Figure  10.-  Continued, 
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(a)  Run  H-17a  of  Reference  9. 

Figure  12.-  Energy  balance  in  the  boundary  layer. 


Drag  parameter,  C,/2,  Re 


Figure  13.-  Composite  plot  of  integral  momentum  equations  results  for 
HIT  data  after  correction  of  blowing  rates  for  energy  balance  results 
v *  0.002. 
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(b)  Det  a  of  Ge*ener  (blowing). 
Figure  18.-  Co v lulled. 


Figure  20.-  Determination  of  the  turbulent  Schmidt  numher  in  the  vicinity  of  the 
wall  from  the  helium  injection  data  of  Reference  17. 


Figure  21.-  Predicted  variation  of  dimensionless  concentration  with  and  without 

transpiration  (Sc  =  0.21,  Sc'  =  0.75). 


Figure  22.-  Predicted  variation  of  dinaniionless  temperature  with  and  without 
transpiration  (Pr  ■  0.7,  Pr '  ■  0.75). 
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Figure  24.-  Comparison  of  the  predicted  dimene lonless  temperature  proiiles  with 
the  data  of  Curl  (Pr  »  0.7,  Pr*  *  0,75). 
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Figure  25.-  Dimensionless  eddy  viscosity  calculated  from  the  data  of 

Kendall  (v  /u  ~  0.0055). 
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(a)  Data  of  Klebanoff,  Reference  SI. 

Figure  26.-  Comparison  of  the  prediction  of  the  entire  turbulent  boundary-layer  velocity 
profiles  with  flat-plate  data,  without  transpiration. 
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(b)  Data  of  Smith  and  Walker,  reference  52. 

Fiqurs  26.-  Continued. 


(a)  Data  of  Kendall  and  Dutton,  Rafarancaa  17  and  16. 

Figure  27.-  Comparison  of  the  prediction  of  the  entire  turbulent  boundary -lays:  velocity 

profile*  with  flat-plate  data. 
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(b)  Oat*  of  Oaaanar,  Rafaranca  11. 
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Figure  28.-  Comparison  of  the  prediction  of  the  entire  turbulent  boundary-layer  temperature 

profiles  with  the  data  of  Curl,  Reference  12. 


(a)  Effect  of  matt  transfer. 

Figure  31.-  Boundary-layer  shape  factor,  H  -  b*/6 
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(b)  Comparison  with  data  of  Reference  50  at  v «  0. 
j  Figure  31.-  Concluded. 
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(a)  Prediction. 

Figure  33.-  Heat-transfer  coefficient  with  and  without  transpiration  (Pr 


Reynolds  number  baaed  on  concentration  thickness,  Pe. 
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Figure  37.-  Effect  of  surface  mass  transfer  on  the  local  skin  friction  at 

a  constant  length  Reynolds  number. 


Figure  38.-  Effect  of  surface  mass  transfer  on  the  local  skin  friction 
constant  length  Reynolds  number  using  modified  blowing  parameter. 
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Figure  40. -  Effect  of  surface  mass  transfer  on  the 
local  mass-transfer  coefficient. 
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Figure  41.-  Effect  of  surface  mass  transfer  on  the  local 

heat-transfer  coefficient. 
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APPENDIX.  A 

THE  ELROD-REICHARDT  CRITERION  FOR  THE  DECAY  OF 
TURBULENCE  NEAP  A  WALL 


Considerable  emphasis  has  been  placed  in  both  the  current  analysis  and 
other  rccent  ane lyses  on  the  continuous  decay  of  turbulence  as  the  wall  xs 
J  approached.  Theoretically,  it  is  possible  to  establish  at  least  one  constraint 
1  with  reqrrd  “o  this  decay.  This  criterion,  herein  referred  to  a-  the  Reichardt 
Elrod  criterion  (Refs.  41  and  47,  respectively),  has  been  derived  for  the  case 
f  of  the  impervious  wall  and  it  is  well  that  it  be  rs -examined  for  the  porous 

wall  case. 

For  the  incompressible  boundary  layer  the  turbulent  shear  is  represented 
j  by  and  iv  is  the  decay  of  this  term  in  the  proximity  of  the  wail  that 

1  is  sought.  This  will  be  accomplished  by  evaluating  successive  derivatives  of 
S  the  term  and  noting  their  values  at  the  wall.  One  supplemental  relation  is 
1  required,  namely,  the  continuity  equation  relating  the  turbulent  components 
1  of  the  three  velocity  vectors 


(A-l) 


for  constant  density  flo^. 

Proceeding  with  the  differentiation  yields  successively 
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The  values  of  these  derivatives  will  depend  on  certain  assumptions  with  regard 
to  the  conditions  at  the  wall.  If  it  is  presumed  that  u'  -  V  -  V  -  0  at 
the  wall,  the  continuity  Equation  (A-l)  requires  that  dv  /dy  *  0.  It  is 
immediately  evident  that  first  and  second  derivatives  of  the  turbulent  shear 
must  vanish,  as  well  as  the  first  three  terms  of  the  expression  for  the  third 
derivative.  The  remaining  term  can  be  rearranged  by  use  of  the  continuity 
equation  as  follows: 
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The  first  term  represents  the  streamwise  variation  of  a  mean  flow  pro¬ 
perty.  At  the  wall  such  variations  are  usually  small  by  comparison  with 
gradients  normal  to  the  wall,  and  certainly  will  vanish  under  truly  one¬ 
dimensional  condition,  for  example,  fully  developed  isothermal  pipe  flow. 

The  second  term  represents  a  correlation  between  derivatives  of  u'  and  w', 
the  two  wall  parallel  components  of  the  velocity.  For  two-dimensional  flow, 
symmetry  arguments  indicate  that  the  correlation  u'w'  must  vanish.  Ifc  is 
invalid,  however,  to  draw  any  conclusions  with  regard  to  the  value  of  the 
correlation  of  the  derivatives.  It  can  only  be  concluded  that  under  conditions 
when  variations  of  mean  properties  with  the  streamwise  coordinate  vanish  at 
the  wall  and  when  two-dimensional  flow  symmetry  exists,  the  third  derivative 
of  the  turbulent  shear  may  vanish. 

It  is  interesting  to  consider  this  result  in  conjunction  with  the  mixture 
length  hypothesis 
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The  first  three  derivatives  of  this  expression  are 
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From  these  expressions  it  is  evident  thit  both  J  and  dl/dy  must  vanish 
at  the  wall  if  either  the  first  two  or  three  derivatives  of  the  turbulent 
shear  stress  are  to  vanish  at  the  wall.  Thus,  the  mixture  length  hypothesis, 
by  its  very  formulation,  eliminates  any  option  with  regard  to  the  value  of 
the  third  derivative  of  the  shear  stress  at  the  wall. 

Therefore,  the  mixing  length  requirements  used  here  are 
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APPENDIX  B 


INTERPRETATION  OF  FRICTION  FACTORS  FROM  PTTOT-PROBE 
READINGS  AT  AND  NEAR  THE  WALL 


It  has  long  been  realized  that  the  experimental  evaluation  of  friction 
factors  under  conditions  of  significant  surface  mass  addition  is  extremely 
complex.  The  momentum  equation 


C,  ,Q  p  v 
f  _  dS  Ko  o 

2  dx  ‘  p1vl 


(B-l) 


although  a  good  tool  when  rates  of  surface  mass  addition  are  low 
(p  vQ/p1v1  <  0.002),  becomes  useless  at  blowing  ratios  in  excess  of  0.003 
for  available  experimental  data.  The  order-of-magnitude  differences  between 
friction  factors  obtained  in  References  9  and  14  at  a  blowing  ratio  of  0.005 
are  ample  deronstration  of  this  fact.  At  this  blowing  ratio,  the  friction 
factor  is  less  than  4  percent  of  the  axial  derivative  of  the  momentum  thick¬ 
ness.  Thus,  three-dimensional  and  time  unsteady  effects  and  errors  associ¬ 
ated  with  the  differentiation  of  experimental  data  becomes  excessive  relative 
to  Cf/2. 

Because  of  this,  an  alternate  technique  was  attempted  in  Reference  17. 
The  technique  involved  the  interpretation  ' f  total  pressure  measurements  in 
a  region  sufficiently  close  to  the  wall  that  the  laminar  sublayer  asserted  c 
significant  influence.  By  evaluation  of  this  influence,  friction  factors 
were  established  hopefully  independent  of  .ny  empirical  postulates  about  the 
character  of  turbulent  momentum  transport. 

The  basic  assumptions  involved  in  thxs  technique  in  Reference  17  were 
two,  r.amely: 

(1)  A  pitot  probe  indicates  a  pressure  equal  to  the  average  of  the 
total  pressures  of  the  undisturbed  incident  flow  over  the  frontal  area  of 
the  probe.  That  is,  if  the  total  pressure  distribution  of  the  incident  flow 
is  represented  by  PT  *  f(y),  then  the  indicated  pitot-probe  pressure  will 
be 


yp+V2 
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PT  dy 
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where  yp  is  the  distance  from  the  wall  to  the  probe  midplane  and  dp  is 
the  probe  thickness. 

(2)  That  the  velocity  (or  total  pressure)  profile  approaches  the  wall 
according  to  the  laminar  sublayer  relation;  which  for  constant  property  flow 
is 
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and  that  turbulence  affects  this  profile ,  at  increasing  distances  from  the 
wall,  in  a  well-behaved  and  reasonably  monatonic  fashion. 

The  validity  of  the  first  assumption  cannot  be  verified  since  no  exact 
calculations  are  known  for  flow  about  a  probe  near  a  wall.  The  principle 
verification  is  based  on  experimental  results  tc  be  presented  in  a  subsequent 
paragraph  of  this  appendix.  The  second  "assumption"  is  logically  acceptable 
without  further  defense  based  on  the  experimental  evaluations  of  the  turbu¬ 
lent  boundary  layers  reported  in  References  48  and  49. 

In  addition  to  the  evaluation  of  these  assumptions,  two  areas  of  evalua¬ 
tion  remain  with  regard  to  the  overall  technique.  These  are: 

(1)  The  accuracies  of  the  input  measurements,  namely,  Pp,  yp,  dp,  and 
v and  their  effects  on  the  interpreted  friction  factor. 

(2)  The  sufficiency  of  the  given  data  and  assumptions*,  that  is,  given 
the  data  and  assumptions,  is  a  unique  value  of  the  friction  factor  clearly 
established? 

It  is  this  final  point  which  is  the  principal  concern  of  the  evaluations 
described  in  this  appendix.  In  Reference  17,  the  establishment  of  the  fric¬ 
tion  factor  was  based  on  hand  fitting  of  the  data  in  light  of  Equations  (B-2) 
and  (B-3)  (see  Figs.  17  through  31  of  that  ref.).  It  has  been  contended 
that  this  process  is  quite  arbitrary  and  thus  does  not  yield  any  valid 
measure  of  friction  factor.  However,  independent  application  of  the  technique 
to  the  same  data  by  different  investigators  has  resulted  in  relatively  con- 
' stent  results . 

In  an  effort  to  further  assess  the  uniqueness  of  the  technique,  several 
analytical  relations  were  postulated  for  the  velocity  (or  total  pressure) 
profile  and  the  arbitrary  constants  (including  C^)  of  these  relations  were 
varied  over  a  fairly  broad  range.  For  each  set  of  constants  the  distribution 

PT  -  f(y> 
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was  obtained  and  the  right-hand  side  of  Equation  (B**2)  was  evaluated  for  the 
probe  locations  corresponding  to  experimentally  measured  values  of  Pp.  The 
sum  of  the  absolute  values  of  the  resultant  errors  was  taken  as  the  measure 
of  merit  of  the  postulated  arbitrary  constants. 

Three  functions  were  postulated.  These  were: 
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Rey  -  a(PT  -  Ps)3  +  b(PT  -  Pfl)2  +  c(PT  -  Pfl)  +  d 
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I  where  P  is  the  static  pressure  and  (P_  -  P  )  *  u*/2  for  constant  property 

5  A  8 

I  flow.1  For  the  first  relation  the  first  derivative  of  the  profile  was  made 
i  continuous  at  the  laminar  limit  Reynolds  number,  two  rela_ 

1  tions  are  based  on  the  Spalding  postulates  (Ref.  53)  for  describing  turbulent 
|  velocity  profiles  and  the  Elrod  constraint  described  in  Appendix  A.2  The  set 
i  of  data  in  Table  B-l  was  selected  (at  random)  from  Reference  17  for  comparison 
with  these  .relations  (Run  C-5*10_3-50,  Station  L) .  The  final  point  was  used 
]  as  an  anchor  point  for  the  postulated  relations,  the  other  four  being  used  to 
evaluate  the  errors  for  a  range  of  values  of  the  remaining  arbitrary  constants. 
For  each  set  of  arbitrary  constants  the  errors  were  determined  at  each  of  the 
four  points,  and  the  sum  of  their  absolute  values  obtained.  This  sum  is 
plotted  in  Figures  B-l  through  B-3  for  the  three  enumerated  relations,  isolat¬ 
ing  the  arbitrarily  assigned  values  of  C^/2  for  the  abscissa.  In  each  case 
a  relatively  distinct  minimum  is  obtained  in  the  range  of  C^/2  between 


1The  u  in  Equations  (B-.  '  ough  (B-7)  should  be  considered  as  equal  to 

u  + yu'u’ ) .  Thus,  tne  constants  of  the  Spalding  type  relations,  for 
example,  would  not  be  expected  to  be  the  same  as  those  reported  in  Ref.  53. 

*  ^Relations  2  and  3  are  consistent  with  the  first  two  and  first  three  deriva¬ 
tives  of  the  turbulent  shear  being  zero  at  the  wall,  respectively. 
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0.00017  and  0.00019.  The  value  reported  in  Reference  17  for  this  station 
was  0.00018. 

Based  on  the  above  results  it  seems  valid  to  accept  the  technique  as 
providing  a  unique  value  of  C^/2  independent  of  the  postulated  profile 
relation.  Therefore,  subsequent  evaluations  will  use  only  one  technique 
(namely,  Eq.  (B-6))  for  the  various  tests  to  be  performed. 

Rather  than  continue  the  rather  awkward  procedure  of  cut  and  try  involved 

in  the  above  technique,  a  least-square  curve  fit  procedure  was  devised  that 

automatically  and  simultaneously  evaluated  those  values  of  a,b  and  Cf/2 

which  produced  the  least-square  error  of  Equation  (B-2).  The  procedure, 

although  iterative,  converged  rapidly  and  produced  excellent  data  fits.  For 

the  case  previously  considered,  Figure  B-4  shows  the  solution  obtained,  its 

integration  according  to  Equation  (B-2),  and  the  comparison  of  this  analytical 

curve  with  the  measured  pitot  pressures. 

« 

A  number  of  tests  were  performed  to  demonstrate  the  sensitivity  of  the 
technique  to  inaccuracies  in  the  requisite  experimental  parameters.  Table 
B-Il  summarizes  these  tests. 

It  is  apparent  from  Table  B-II  that  no  major  amplification  of  the  errors 
in  the  basic  measurements  -'ccurs  when  the  friction  factor  is  evaluated.  This 
is  in  contrast  with  the  momentum  equation  technique  where  violent  error 
amplification  occurs  at  high  rates  of  surface  mass  addition. 

Having  adopted  an  analytical  expression  for  the  velocity  profile  the 
next  step  in  this  evaluation  study  was  to  reevaluate  (with  the  least- square 
curve  fit  program)  all  the  friction  factors  reported  in  Reference  17.  These 
values  are  presented  in  Tab]'*  "-III.  These  values  differ  only  slightly  from 
the  originally  reported  valuer. 

Probably  the  most  significant  test  of  the  overall  technique  is  provided 
by  comparison  of  values  of  C^/2  obtained  without  mass  addition  with  those 
generally  reported  in  the  literature  and  those  obtained  by  use  of  the  momen¬ 
tum  equation  using  the  data  of  Reference  17.  This  comparison  is  shown  in 
Figure  16  where  it  can  be  concluded  that  the  procedure  does  indeed  perform 
a  creditable  job.  It  is  significant  to  note  that  with  increasing  surface 
mass  addition,  the  total  pressure  gradient  is  significantly  reduced  near 
the  wall  and  thus  the  validity  of  Equation  (B-2)  becomes  less  susceptible  to 
question.  That  is,  if  the  technique,  and  by  inference  Equation  (B-2),  is 
valid  when  there  is  no  surface  mass  addition,  it  should  still  be  valid  when 
there  is  surface  mass  addition. 


B-5 

It  would  seem,  based  on  the  evaluations  described  above,  that  friction 
factors  can  be  obtaired  through  interpretation  of  pitot- probe  measurements  at 
and  near  the  wall,  and  further,  that  the  technique  involved  does  not  degener¬ 
ate  as  rates  of  surface  r  ass  addition  are  increased. 
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TABLE  B-I 


DATA  USED  TO  TEST  FRICTION  EVALUATION  PROCEDURES 
UTILIZING  PROFILES  (B-4)  THROUGH  (B-7) 


Rey 

BW 

ill 

197 

0.0086 

197 

0.0086 

270 

.0172 

393 

.0346 

516 

.0493 

761 

.0702 

TABLE  fl-II 


SENSITIVITY  OF  FRICTION  FACTOR  TO  INACCURACIES 
OF  PROFII£  PARAMETER 


Nominal 

Perturbation 

Friction  Factor 

r actor 

Value 

Value 

Chanqe 

(%f 

Value 

Change* 

(%) 

V 

</»» 

0.00538 

0.006 

11.5 

0.000168 

-  9.2 

Probe  height 

.16 

.018 

12.5 

.000159 

-14.1 

Scale  shift  in  y 

0 

.002 

— 

.000131 

-29.2 

0)  0) 
u  c 

0.008  in. 

.007 

.009 

28.6 

.000212 

14.6 

3  (0 
(Q  4J  | 
t0  Q, 

.011  in. 

.014 

.016 

14.3 

.000207 

11.9 

<u  <u  0( 
u  N 
a 

.016  in. 

.028 

.030 

7.1 

.000189 

2.2 

<0  C  V 

A  <o 

0  • 

.021  in. 

.040 

.042 

5.0 

.000171 

-  7.6 

u  c 

.031  in. 

.057 

.059 

3.5 

.000190 

2.7 

♦Nominal  value  of  friction  factor  ■  0.00185. 


TABLE  B-III 


FRICTION  FACTORS  INTERPRETED  FROM  PITOT  PRESSURE  PROFILES 

ADJACENT  TO  THE  WALL 


Run 

Station 

Cf/2 

(X103) 

Vui 

(X103) 

Ree 

Rex 

(xl0‘y) 

From 

Ref.  17 

Least- 

Square 

Fit 

0 

2.5 

2.77 

0.1098 

686 

0.130 

E 

2.4 

2.81 

.1192 

931 

.249 

G 

1.75 

1.76 

.1076 

1332 

.490 

H 

1.95 

1.89 

.1133 

1730 

.634 

C-iu  -bu 

I 

1.6 

1.67 

.1140 

2011 

.893 

J 

1.75 

1.75 

.1203 

2433 

1.083 

K 

1.45 

1.56 

.1203 

2820 

1.375 

L 

1.6 

1.58 

.1184 

3377 

1.658 

E 

2.35 

1.129 

1063 

.248 

G 

1.49 

1.122 

1573 

.490 

I 

HdjH 

1.25 

1.129 

2416 

.887 

I* 

1.30 

1.120 

2472 

.895 

C  — 1U  -DU 

K 

1.17 

1.118 

3581 

1.388 

L 

1.51 

1.032 

4308 

1.678 

M 

BSSH 

1.13 

1.080 

43?  ^ 

2  .067 

N 

nfl 

1.08 

1.076 

5839 

2.290 

E 

1.47 

i  .97 

2.22 

1237 

.252 

G 

1.21 

1.18 

2.17 

1903 

.499 

H 

1.1 

1.30 

- .  16 

2533 

.695 

I 

1.05 

.92 

2.19 

2990 

.906 

C-2Xl0”3-50 

J 

.96 

1.10 

2.15 

3704 

.102 

K 

.91 

.82 

2.30 

4350 

i  312 

L 

.89 

.94 

2.13 

5435 

1.670 

M 

.84 

.80 

2.19 

6314 

2.002 

N 

.80 

.82 

2.21 

7482 

2.304 

E 

1.25 

1.55/1.61 

3.23 

.248 

G 

1.0 

1.18 

3.25 

.492 

H 

.90 

.98 

3.21 

.681 

J 

.64 

.65 

3.27 

4510 

1.083 

C-3X10  J-50 

K 

.58 

.54 

3.26 

5453 

1.381 

L 

.58 

.63 

3.26 

6725 

1.665 

M 

.46 

.44 

32.5 

7886 

1.987 

N 

.46 

.48 

3.22 

9506 

2.230 

E 

.68 

0.86/*.  90 

5.43 

1740 

.254 

G 

.41 

.49 

5.42 

2960 

.504 

H 

.36 

.42 

5.40 

4020 

.698 

C-5xl0“3-50 

J 

.21 

.19 

5.57 

6170 

1.110 

K 

.18 

.19 

5.42 

7740 

1.410 

L 

.18 

.19 

5.38 

9800 

1.700 

M 

.16 

.18 

5.37 

13S0C 

2.037 

Range  of  values  resulting  from: 


0.69  <  2-  <  0.71 


45/  -i-  ■■  ^  y  7  10‘S 

*  <  qx  dRe  <  7-10 

*  J 


at  Re  750 

y 


with  optimum  Re1<uainar  limit 


Friction  factor,  C^/2,  X104 


Figure  B-l.-  Effect  of  varying  arbitrary  constants  of 
Equation  (B-5)  on  the  resultant  error  of  Equation  (B-2) 
for  the  measured  data  points  of  Table  B-l. 


Sum  of  absolute  values  of  error:  ,  xlO 


1  dq 


Curve  obtained  by  setting 
-3-  “  0.69  ^ 


4.5  < 


_1_  dq 

qi  dKev 


<  7X10"5 


at  Rey  -  750 


with  Cf/2  optimized  at  each 
value  of  (l/q1) (dq/dRey) 


Friction  factor,  C^/2,  XlO4 


Figure  B-2.-  Effect  of  varying  arbitrary  constants  of 
Equation  (B-6)  on  the  resultant  error  of  Equation  (B-2) 
for  the  measured  data  points  of  Table  B-l. 
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Figure  B-3.-  Effect  of  varying  arbitrary  constants  of 
Equation  (B-7)  on  the  resultant  error  of  Equation  (B-2) 
for  the  measured  data  points  of  Table  B-l. 


Figure  B-4.-  Technique  for  establishing  local  skin  friction  from  dynamic  pressure  measurements 

throuoh  least  square  fit. 


